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ABSTRACT 


A second order numerical method employing reference plane 
characteristics has been developed for the calculation of 
geometrical ly complex three dimensional nozzle-exhaust flow 
fields, heretofore uncalculable by existing methods. The 
nozzles may have irregular cross sections with swept throats 
and may be stacked in modules using the vehicle undersurface 
for additional expansion. The nozzles may have highly non- 
uniform entrance conditions, the medium considered being an 
equilibrium hydrogen-air mixture. The program calculates 
and carries along the underexpans i on shock and contact as dis- 
crete discontinuity surfaces, for a nonuniform vehicle ex- 
ternal flow. Additionally, shock formation due to coalescence 
is detected. A wide variety of geometric problems may be 
considered since the reference plane method has been developed 
for three separate coordinate systems, all incorporated into 
a single program. 
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I . INTRODUCTION 

The propulsion system of a scramjet powered hypersonic aircraft 
is one of the dominant portions of the vehicle structure, and 
renuires a high degree of veh icl e-propul s ion system integration 
to achieve optimum performance. The engine exhaust flow, be- 
cause of physical area limitations, may be underexpanded at the 
nozzle exit and, in order to obtain maximum propulsive efficien- 
cy, the vehicle afterbody under-surface is used to provide addi- 
tional expansion. This results in a three dimensional nozzle 

r 

flow whose boundaries are defined by the solid walls of the noz- 
zle and the vehicle surface and by the stream surface separating 
the nozzle flow from the external flow field of the vehicle. 

The effect of this complex interaction of the nozzle flow with 
the external flow, in terms of the forces and moments produced 
on the vehicle, can greatly influence the vehicle stability and 
aerodynamic performance, as well as the engine performance and 
must be considered in design studies of hypersonic airbreathing 
aircraft. 

While much effort has been expended on developing numerical tech 
niques for the solution of three dimensional supersonic flow 
fields, only a few working computer programs have evolved from 
these efforts. The numerical schemes generally fall into one 
of the following categories: (1) finite difference methods; 

(2) bi -character istic methods; and (3) reference plane methods. 
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The reference plane method has been selected for this effort 
due to the geometric complexity of the flow field and the oc- 
curence of both embedded and external shock waves. The finite 
difference approach becomes difficult to apoly when the bound- 
aries cannot be reoresented by simple analytic functions and 
the accuracy of shock calculations tends to be very poor if the 
shocks are not steeply inclined with respect to the marching 
direction. A bi-characteristic approach may lead to problems 
of complex grid ordering and interpolation procedures and the 
calculation of internal discontinuities can be a formidable 
task. The reference plane method can alleviate these difficul- 
ties provided that a coordinate system is chosen which is suit- 
ably aligned with the flow field. Additionally, in the local 
reference planes, there is a strong analogy to a two dimensional 
flow field, hence the logic and numerical grid ordering appro- 
priate to a two dimensional system can be applied. This is 
particularly important for the geometrically complex flow fields 
being considered, since the reference planes can be locally 
oriented to conform to the vehicle geometry, thus leading to ex- 
cellent flow visualization. 

Under this contract two families of exhaust nozzles will be an- 
alyzed. The first family (Vehicle I) will be characterized as 
having nozzle modules with cross-sections which are rectangular 
in shape. These nozzles nay be arranged in multiples and dis- 
charge into a common nozzle, as is the case, for example, of 
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several engine exhaust nozzles using a vehicle undersurface 
for additional expansion area. This arrangement is shown in 
Figure (1). The second family (Vehicle II) is characterized 
as having a throat of curvilinear shape with one dimension 
much larger than the other. The throat may be sweptback in 
the downstream direction. This configuration is depicted in 
Figure (2). 

The reference plane method was first proposed by Ferrari in 
1949 (Reference 1) and further discussed by Ferri in (1954) 
(Reference 2). More recently, the method was used by Moretti 
(References 3 and 4) and Rakich (References 5, 6 and 7) for 
the analysis of external flow fields about inclined bodies of 
revolution. 

The choice of the reference plane system to be employed for 
external flow fields around inclined bodies of revolution is 
quite clear, namely an r, 6, z system where z coincides with 
the body or wind axis and the cross flow is in the 9 direction. 
For the analysis of scramjet exhaust nozzle flow fields the 
choice depends on the specific nozzle geometry; different con- 
figurations require different reference plane systems. 

In addition, for a given configuration different reference sys- 
tems may be used in different regions of the flow field so that 



v_ ENGINE 
MODULE 



FIGURE 1. VEHICLE I, CONFIGURATION 
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each system is locally aligned with the flow field and the 
boundary conditions can be properly satisfied. 

In the reference plane method, the configuration of the refer- 
ence planes is chosen such that the primary flow variations 
occur within the reference planes while the coordinate normal 
to them describes the local cross flow. The governing equa- 
tions are written in this reference coordinate system with the 
terms representing derivatives normal to the reference planes 
treated as forcing functions. The characteristic directions 
are those of the two dimensional system in the reference plane 
while the compatibility relations differ from those of the two 
dimensional system by the inclusion of the cross derivative 
terms, which are evaluated on the initial surface at each 
marching step. Note that flow variations must be continuous 
with respect to the cross flow direction, while wave-like dis- 
continuous flow is permissible in the reference planes. 

The gas mixture considered consists of air and combustion prod- 
ucts (with hydrogen as the fuel) and is assumed to be inviscid 
and in chemical equilibrium. Three parameter curve fits (p, $, 
h) have been developed from the data of Reference (8) for the 
necessary thermodynamic properties, as described in Appendix I. 
The flow may be rotational and nonhomentropi c . 

The numerical procedure in a given reference plane employs a 
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Hartree type numerical grid* the streamline projections being 
traced in the individual reference planes. This approach af- 
fords excellent flow visualization through the tracing of these 
quasi streamlines, and provides a self enforced mesh tightening 
in regions of flow compression, since the quasi streamlines 
converge in these areas, thereby decreasing the axial step size. 
This also provides an aid in the detection of embedded shock 
waves . 

The overall numerical approach is second order both with respect 
to reference plane characteristic calculations and the evalua- 
tion of derivatives. That is difference coefficients appearing 
in the compatabi 1 i ty relations as well as derivative functions 
are averaged during a global second order computation. The 
program, however, does provide for a purely first order computa- 
tion as well as local second order characteristic calculations 
using only the initial value derivatives. During global cal- 
culations shock and contact surface cross angles are updated 
in addition to cross derivatives. 

All interpolations are linear since this always led to stable 
numerical results as well as ensuring convergence of the dif- 
ference solution to the exact solution. However, care must be 
exercised in using an extremely course grid spacing in the re- 
ference planes since excessive numerical diffusion of the solu- 
tion will resul t . 
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1 1 . GOVERNING EQUATIONS 


The equations governing the steady three-dimensional flow of 
an inviscid gas mixture in chemical equilibrium are as follows: 

Continuity 


V • (pv) = 0 


( 1 ) 


Momentum 

p(7*V)7 + 7p = 0 (2) 

Conservation of Stagnation Enthalpy 

7 ♦ V(H) =0 (3) 

Conservation of Entropy 

7 ♦ Vs = 0 (4) 

Equation of State 

(!&)*!£. a 2 (5) 

' Sp ' s p 

Constancy of Equivalence Ratio Along Streamline 

7 • 7$ = 0 ( 6 ) 

Curve Fit For Isentropic Exponent 

r = f ( h , p , $ ) * 


♦ Fits Tor r have beerf expressed in polynomial form from properti es 
tabulated in Reference (8) as described in Appendix I. 
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Combining Equations (4) and (5) yields the relation 
v . Vp = a^Vp 

and hence the continuity equation may be written in the form 

7 . Vp + a 2 pV .7=0 (8) 

The analysis of the nozzle-exhaust flow fields for Vehicles 
I and II has required the use of three types of coordinate sys- 
tems. 

(1) Line Source System: This is an r-9-z system 

with the flow predominantly aligned with the r 
direction as depicted in Figure (3). The refer- 
ence planes are the planes 6=constant and the 
marching direction is r. The velocity vector 

in this system is depicted in Figure (4). 

(2) Cartesian System: This is an x-y-z system with 

the flow predominantly aligned with the x direc- 
tion as depicted in Figure (5). The reference 
planes are the planes y-constant and the march- 
ing direction is x. The velocity vector in 
this system is depicted in Figure (6). 
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FIGURE 4. 


VELOCITY VECTOR IN LINE SOURCE SYSTEM 
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AAA 

V«Ul x + Vly + Wl 2 
() = V cos i|) Q = ( u^+w^ ) 2 

u = Q cos <f> = tan -1 {w/u) 

v = Q tan i|> \p = tan - *(v/Q) 

w = Q sin 


FIGURE 6. VELOCITY VECTOR IN CARTESIAN SYSTEM 
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(3) Cylindrical System: This is an r-6-x system 

with the flow predominantly aligned with the 
x direction as depicted in Figure (7). 

The reference planes are the planes 6=constant and the march- 
ing direction is x. The velocity vector in this system is 
depicted in Figure (8). 

In the overall nozzle-exhaust analysis, let us define the in- 
ternal flow as that portion of the flow field bounded by the 
vehicle undersurface, cowl and sidewalls (or planes of sym- 
metry), and the external flow as that portion of the flow field 
bounded by the vehicle undersurface and an outer prescribed stream 
surface (or planes of symmetry). 

In analyzing the internal flow field for Vehicle I, a line 
source system would be employed if the sidewalls were curved or 
aligned with a ray other than one parallel to the x direction 
as shown in Figure (3), while a Cartesian system would be employed 
if the sidewalls are parallel to the x-z plane, as in Figure (5). 
For the internal flow portion of Vehicle II, a cylindrical 
system is employed as shown in Figure (7). 

If a central module of Vehicle I is being analyzed, the side- 
walls are replaced by planes of symmetry in the external flow 
and the Cartesian system shown in Figure (9a) is employed while 




FIGURE 7. CYLINDRICAL SYSTEM 


V 



v i e + w i r 

h 

Q = (u 2 +w 2 ) 2 
= tan _ 1 (w/u) 
if) = tan"*(v/Q) 


II 


I CAL SYSTEM 



FIGURE 9a. EXTERNAL FLOW -VEHICLE I CENTRAL MODULE 
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for Vehicle II the external flow is computed as indicated 
in Figure (9b). The end modules are analyzed only for Vehicle 
I and are computed separately in a locally Cartesian system as 
depicted in Figure (10). 


A. Line Source System In this system, the refer- 
ence planes are the planes 0 = constant and variations in 
these planes are expressed by differentiation in the r and 
z directions. The primary criterion in selecting a refer- 
ence plane system is that the reference planes be aligned 
such that the velocity component normal to them is as small 
as possible in order to obtain the maximum axial step size 
in the difference scheme. 

The scalar equations of motion are: 


Continuity 


pa 2 (u + -£ + — + w ) 


r r r "z' + up r + ~r^ + W P Z 


v Pe + 


= o 


r - momentum 


p(uu, 


VUq 

— + — - + wu_) + p„ = 0 
Y » y z ' r 


(9) 


( 10 ) 
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EXTERNAL SIDEWALL 
CALCULATION IN Y t Z SYSTEM 


FIGURE 10. EXTERNAL FLOW END MODULE 
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6 -momentum 


p(v^. + W v z + uv r + ^) + p- ■ 0 

z-momentum 

, VWQ 4 

P ( ww z + uw r + — ) + p z = 0 
The following relations are valid along streamlines 

d£ _ r !P 
P p 

H - constant 

$ » constant 


(ID 


( 12 ) 


(13) 


(14) 

(15) 


a 2 = IP (16) 

P 

The system of Equations (9) - (16) constitutes a system of 
quasi-linear partial differential equations in the variables 
u, v, w, p, p, H, $ and a. The system can be solved 
numerically by an ordinary finite difference formulation based 
on arbitrary coordinate directions or by a finite difference 
scheme based on characteristic directions. The ordinary finite 
difference formulation leads to a great simplification of the 
problem if the flow is everywhere continuous and the boundaries 
are smooth. However, major problems in accuracy are encountered 
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in devising means of satisfying boundary conditions on surfaces 
not corresponding to coordinate planes, which then necessitates 
interpolation or extrapolation. Also, even on boundaries cor- 
responding to coordinate planes, the derivatives at the boundaries 
can only be approximated by one-sided differences. 

A true three-dimensional characteristic calculation is quite 
complex from a computational point of view. The usefulness of 
the characteristic method is weakened when three independent 
variables are considered, since the system of equations no 
longer reduces to a system of ordinary differential equations 
but rather to a system of partial differential equations in two 
independent variables on characteristic surfaces. The com- 
plexities involved in analyzing nonuniform flow fields with 
embedded shocks and complicated geometries, makes the use of a 
three-dimensional characteristic approach numerically complex 
and far less attractive than the present method. 

The method employed to numerically solve the system of 
Equations (9) - (16) in this analysis uses a finite difference 
grid consisting of the characteristic directions obtained by 
considering the momentum equations in the r and z directions 
and the continuity equation. These characteristics are not 
the bi character! sti cs of the three-dimensional system, but rather 
represent the projection of the bi -character! sti cs onto the refer- 
ence planes. 



TR 166 - Vol. I 


Page 22 


The r-nomentum, z-momentum and modified continuity equation are 
rewritten by cutting the terms involving derivatives in the 9 
direction onto the righthand side. These ecuations now take 
the form 


p(uu r + wu z ) + P r = 


p(uw r + WW Z ) + P z 




(17) 


(18) 


pa 2 (u r + w z ) + up r + wp z 


= -a 


» /U 

■P ( r 


— ) - 
r ' 


v Pe 

r 


(19) 


The initial value plane is one of constant r in this coordinate 
system and the righthand side terns, including the cross deriva- 
tives, will be evaluated in this initial nlane. Then, the system 
(17) - (19) is a system of 3 ouasi-linear partial differential 
equations in the variables u, w and p and the lefthand side is 
equivalent to the corresponding two-dimensional system in the 
r-z plane. The characteristic directions for this system are 
those of the two-dimensional system and the compatibility rela- 
tions differ from the two-dimensional relations by the inclusion 
of cross derivative terms in the 9 direction. 

Letting, a, $ and ^ represent the velocity vector in nlace of 
u, v and w (as shown in Figure 4)» the velocity vector v 
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may be wri tten , 


A A A 

v = u i r + v ig + w i f = 

( 20 ) 

A A A 

q (cos<t> i r + tan ^ ig + sin<j> 1 Z ) 


where q represents the projection of the velocity v onto the 
reference plane 6=constant. 


The characteristic directions are given by the equations 


= (q/a) 2 cos £ sin Jut-1 = Xi (21) 

, 'j 2 

c - (q/a) 2 cos 2 <p - 1 

where $ 2 = q 2 /a 2 - 1 
and 


(dz) = tan <P s w (22) 

S.L. u 


which represents the streamlines projection onto the reference 
plane . 
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Along a C± characteristic, the compatibility relation takes 
the form 

pq 2 d <t> ± 0 dp = ( F ± ) dr 

where 

r 

F , = ! ( s i n <j> 

± r I 

I 

U 

+ liL-i p 0 + cos (J> ] - tan p (4» 0 ) 
pa 2 i 

(sin X + + cos d> ) - X ± tan 2 tpl 


- cos <p Xj 


I 


i 2 
\ COS 


( 23 ) 


(24) 


B. Cartesian System In this system, the reference 
planes are the planes y-constant and variations in these planes 
are expressed by differentiation in the x and z directions. 
Analogous to Equations (9) - (12) for the line source system, 
the scalar equations in the Cartesian system take the form: 


Conti nui ty 

pa 2 (u x +v y +w z ) + up x +vp y +wp z =0 (25) 

X-MlOfflSntlliT! 

p (uu + V U +WU ) + p =0 
' x y z x 


( 26 ) 
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y-momentum 

p(uv x +vv y +wv 2 ) + p y =0 (27) 

z-momentum 

p{uw x +vw y +ww 2 ) + p z =0 (28) 

The streamline relations (Equations 13 - 16) are valid for all 
coordinate systems. The x-momentum, z-momentum and continuity 
equations are rewritten by putting the terms involving deriva- 
tives in' the v direction onto the right hand side. The initial 
value plane is one of constant x and the terms on the right 
hand side are evaluated on this plane. 


Letting q, tp and represent the velocity vector as depicted in 
Figure (6), the characteristic directions are given by the equa- 
tions 


where 


dz 

dx 

C± 


(q/a) cos p sin d> ft _ .± 

P ? A 
(q/a) cos <p - 1 


e 2 = q 2 /a 2 


1 


(29) 


and 


dz 

dx 

S.L 


tan <j> = H 
u 


( 30 ) 


Along a C+ characteristic, the compatibility relation takes the 
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form 


pq 2 dt ± 6 dp = (F + ) dx 


(31) 


where 


F + = po' 


(sin <j> - X* cos <J> ) 


v y tan ^ 

— ? + — ? M 

cos^ pa " 1 


- tan \J>(4> v ) (cos <J> + X + sin <j> ) 

V 


(32) 


C. Cylindrical System In this system, the reference 
planes are the planes ©^constant and variations in these planes 
are expressed by differentiation in the x and r directions. Anal- 
ogous to Equations (9) - (12) for the line source system, the 
scalar Equations in the cylindrical system take the form: 

Conti nu i ty 


pa^(u + — + w + — ) up + ~p Q +wp 

M v x r r r' p x r ^0 y r 


= 0 


(33) 


x-momentum 


p(uu x + ^ u e + wu r ) + p x = 0 


(34) 


©-momentum 


p(uv x + £ V, + 


+ uv r ) + ^ = 0 


wv 

r 


(35) 
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r-momentum 

2 

p(uw r + — Wg " + ww 2 ) + p r - 0 (36) 

The x-momentum, r-momentum and continuity equations may be re- 
written with the terms involving derivatives in the 0 direction 
being put on the right hand side. Letting q, <j> and represent 
the velocity vector as depicted in Figure (8), the characteristic 
directions are given by the equations 


dr 

dx 


(q/a) 2 sin cos <f> 
2 2 

(q/a) cos - 


(37) 


where 



and 


dr 

dx 


S.L. 


w 

= tan <j> = - 


(38) 


Along a C + characteristic, the compatibility relation takes the 
form 


pq^ d<f> ± $ dp = F dx 


(39) 


where 


F 1 = £ 


q{sin <J> - X cos <f>) 


v e 


vq e 


( 40 ) 


+ If® + q sin <j> j - qv <j> B (A* sin <J) + cos <}>) + v2 
pa 2 
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III. NUMERICAL SCHEME 

Consider the numerical grids depicted in Figures (12a), (12b) 
and (12c) for the line-source, Cartesian and cylindrical sys- 
tems respectively. The reference planes are the planes J= 
constant and in a given reference plane, the index I goes from 
1 (lower wall) to IMAX(J) (cowl or boundary stream surface). 
The reference plane 0=1 is always a plane of symmetry. 

The variables p, q, <t> , ip , p, h and $ are specified at all the 
mesh points I, 0 in the initial surface (r=rj for line-source 
system in Figure (12a); x=Xj for both the Cartesian and cylin- 
drical systems in Figures (12b) and (12c) respectively). It 
is desired to calculate properties at the station (for the 
line-source system) or (for the Cartesian and cylindrical 
systems) where the step-size Ar or Ax is constrained by stabil 
ity considerations to be discussed in a subsecuent section. 
Note that a uniform step-size is taken for all grid points 
since the evaluation of cross derivatives on the surface r 
or x-constant is required. 

A. Cross Derivatives at Interior Points Before pro- 
ceeding with the calculation, cross derivatives must be evalu- 
ated and stored at all mesh points I, J on the initial surface 
This is complicated by the fact that the spacing of the grid 
points is arbitrary in each reference plane J and that the 
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* t 





TR 166 - Vol . I 


Page 33 


spacing (as or Ay) of the reference planes is arbitrary. It 
is required that each reference plane contain an identical 
number of points. Referring to Figure (11) all cross deriva- 
tives (af/ay) are computed by combining derivatives (af/az) 
and derivatives (af/as). This requires that there be a one to 
one correspondence between points (1,0-1), (I,J) and (I,J+1) 
which form the arbitrary surface n = constant. This approach 
has the advantage of eliminating unnecessary interpolation of 
initial data. In addition, derivatives at discontinuities, 
fixed boundaries and interior points are computed in a consis- 
tent manner. 

The relation that exists between derivatives in the z,s and 
y(or e) directions on the initial surface is: 



where as is defined as 

as = [(Ay ) 2 + (az) 2 ] 2 

Thus, knowing the derivatives (af/as), and (af/az), as well as 
the direction angles (ay/as) and (az/as) allows us to compute 
the cross derivative (af/ay) at every point (I,J) of flow field 


(41) 
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(1) Line Source System: At all interior mesh 

points (1,0) the derivative ( a f / a e ) is required for the 
property f. Referring to Figure (12a) the n = constant sur- 
face is formed by points (I,J-1), (I,J) and (I,J+1). Then 


AS i " ^ 2 i,j " 2 1 ,0- 1 ^ 2 + r2 ( 0 j ' e j-i^ 2] 


(42) 


A5 2 " [ ^ Z I,J+1 ' Z I,J^ 2 + r2 ^ 0 J+l * 0 J^ 2] 


using the standard centered difference formula for nonuniform 
spaci ng 


36 

e J+1 

AS, 

i AS, A S o 

- e J ( — - - — -) 

V AS 2 ASj ; 

- 9 J_1 

AS 0 

(.s') 

3 S 

J + i 

<n*> 

(ASj + AS 2 ) 

AS, AS~ 

> Z J ( L £) 

I 'AS 2 ASj' 

- z J_1 

AS, 

A! 

Z I 

Z I 

3 S 



( as x + A S g ) 




and for any property f 

f J+1 (- i) 

3 f _ 1 (as 2 J 

. f j . th) . fj-i ( as 2) 

T I ^AS 2 AS j r I 'ASj' 

3 S 

(ASj + as 2 ) 


(43) 


( 44 ) 


AZ 1 = Z I,J " Z I-1,J 
Az 2 = Z I+1,J " Z I,J 


using 
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then 


ll 

dz 


1 Az i 

f J ( t) 

T I+1 'hZ z } 


, AZ, AZ-3 

f J ( L L\ 

T I k AZ 2 AZ ^ 


1 AZ ? 
f J { L\ 

r i-l AZ j ' 


(AZj + AZ 2 ) 


( 45 ) 


the cross derivative can thus be determined from Equation (41) 


/ 5f v _ r 9_f /9f\ 9z_ , ,96 

^ 96 ' “ L 9S " 9 Z 9 S J/ 9S 


(46) 


2. Cartesian System: The derivative (9f/9y) for all 

properties f are required. The procedure is analogous to that 
described for the line source system except that referring to 
Figure (12b) 

AS 1 = " z I , J - 1 ^ 2 + ^ y J " y J-l^ 1 

(47) 

2 2 h 

As 2 = [ ^ Z I,J + 1 ” Z I , J ^ + ^ y J + l ' y J^ 1 
The cross derivative is then determined from Equation (41) 


(il) = 

' 9y ^ 


[ — 

9 S 


( 11 ) 

V 9 Z 7 


9 Z 
9 S 


]/iX 

9 S 


( 48 ) 


3. Cylindrical System: The derivative (9f/90) 

is required at all interior mesh points (1,0) for the proper- 
ties f. The procedure is analogous to the line source system 
except that referring to Figure (12c) 
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AS 1 = [{r I,J ' r I,J-l )2 + ?2 ^ e J • e J-l )2]2 


as 2 s C ^ r I,J + l " r I,p 2 + r2 ^ 9 J + l * 


(49) 


and the cross derivative is from Equation (41) 


, 3f v _ . 3jf /3f\ 3_r . ,3J0 

1 as l 3r' as J/ as 


(50) 


B. Cross Derivatives At Boundary Points - The cal- 
culation of cross derivatives at discontinuity surfaces and 
walls are formally similar to the calculation of cross deriva- 
tives at interior points of the flow field. The arbitrary n = 
constant surfaces are replaced by physical walls, shock or con- 
tact surfaces. However, care must be exercised in the com- 
putation of ( a f / s z ) derivatives. Due to stability considera- 
tions, forward and backward differencing for these derivatives 
replaces the usual centered differences. 

As shown in Figure (13) no distinction is made between an upper 
wall boundary or the lower surface of a shock or contact dis- 
continuity. Similarly, no distinction is made between a lower 
wall boundary and the upper surface of a shock or contact dis- 
continuity. 

For example, in a Carteisan system (9f/3y) the required cross 
derivative at the boundary is defined by 
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a 

I - Y 




FIGURE 13. UPPER BOUNDARY OR LOWER SURFACE 
OF DISCONTINUITY 



TR 166 - Vol . I 


Page 38 


il - r if. /if\ 1L i ill 

ay 1 as ^az' as J/ as 


(51) 


where (af/az) is replaced by the appropriate forward or back- 
ward difference formula. 


Referring to Figure (13) a combined differencing scheme for 
(af/az) on either a lower type boundary (L=l) or an upper type 
boundary (l=-l) is 


& 2 f I+2L~ (1+a)2 f I+l t (lt2a) f l 
AZj (1+A) 


(52) 


where 

42 l = 2 I - Z I+L 

and 

A = (AZj)/(Zj + l - 2 1 + 2|_ ) 


C. Step Size Criterion - The step size is limited 
by the Courant-Fri edri chs-Lewy (C-F-L) stability criterion 
which states in geometric terms that the domain of dependence 
of the differential equations must be embedded within the domain 
of dependence of the difference scheme. Referring to Figure 
(14) for the (1,0) point depicted, the domain of the difference 
scheme is enclosed by the points (I,J-1), (1+1,0), (I,J+1) and 
( 1 - 1 , 0 ). 
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Then, the C-F-L stability criterion states that the region formed 
by the intersection of the Mach fore cone through I with the plane 
x=x i (the shaded ellipsoidal area in Figure 14) must be enclosed 
by the domain of the difference scheme. This criterion then allows 
one to determine the maximum marching step (aXj j = x i j “ x i j) 
for each mesh point I,J and hence select the overall marching step 
ax^g as the minimum of all the AXj j. 

When discontinuities are present in the flow, the step-size is re- 
stricted by the consideration that no derivatives may be made by 
crossing a discontinuity. Referring to Figure (15) derivatives 
in the z direction at point I cannot be made by the usual centered 
difference formula since this would entail crossing a discontinuity 
and hence violate the C-F-L stability criteria. Thus, as previous- 
ly described in Section B above, the appropriate forward or back- 
ward difference is used in making the z derivative. The domain of 
dependence in this case is given by the ellipsoid of shaded area 
of Figure (15). Thus, the C-F-L condition at a boundary is ap- 
plied in a formally analogous fashion to that of an interior 
point. 


D. Interior Point Calculation* - Referring to 
Figure (11), a typical interior point calculation for the line 
source system in the reference plane J is performed as follows. 

*Let C refer to the point ( I , J ) and D refer to the point (1,0) 
in Figures (11), (12) and (13). 



Page 41 


DISCONTINUITY 

SURFACE 


o GRID POINTS 


MAIN OF 
FACE 


DEPENDENCE A+A 


TR 166 - Vol. I 


Page 42 


The point C is located by projecting a quasi -streaml i ne from D 
a distance Ar downstream (Ar being the overall marching step). 
Using Equation (22) in difference form 

s z D + tan <j> D Ar 

Points A and B are located by shooting back characteristic pro- 
jections from point C, locating these points between the mesh 
points 1-1 and I, and I and 1+1 respectively, by an iteration 
procedure, using Equation ( 21 ) in difference form: 

2 A * 2 C - U) A 4r 
B B 

All required properties can be obtained at A and B by a linear 
interpolation, including cross derivatives which are stored at 
the grid point locations. 


(53) 


(54) 


The pressure and flow deflection angle in the reference plane 
<f >0 can be determined by solving the linear relations along AC 
and BC. 


(pq 2 ) A (* C - *A> * S A < P C ' P A> = F A lr 
*„) - S B (p c - P B ) ■ F' Ar 


<pq 2 ) 


B *B' R B Vh X " K B' ” ' B 

The 9-momentum equation (Equation 11) may be rewritten in the form 


p u { v + — v ) = 
v r u z ' 


^9 pvv e P uv 


si nee 


ar az H as 


(55) 

(56) 

(57) 

(58) 
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(where — represents the change along the direction CD, a 

t) s 

streamline projection), we obtain the difference eouation: 


V C s V 1 


r cos <ft ^ 


Ha 

- — + tan ili tan \p (q6 ) 

po \ 


+ + q cos 4> 

cos^ 


D 


Ar 


All streamline relations are valid along the actual stream- 
line DC. Properties at D* may be obtained to first order by 
the relation 

f lT* = f D “ ^ f 9^D A0 

which upon incorporating streamline geometry reduces to 


(59) 


(60) 



* f fpD 


',, tan t 1 Ar 
r cos d> j 


The pressure density relation, expressed in differential form 
along a streamline in Equation (13) may be integrated holding 
F constant for the small integration step CD*, yielding 


(61) 



( 62 ) 
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The Bernoulli relation 


dp dV L 

+ -5- = 0 
P 2 


be integrated along D*C, making use of Equation (56), 


yielding 


v = 
C 


vL - 2i_ 

v d* r-i 


The conservation of total enthalny states that h^-h^* and 


hence , 


h c ' h D* + X Cv*. - V?) 


2 D* C 


For inviscid flows, the equivalence ratio remains constant 
along a streamline, hence, 


C D* 


The i sentropic exponent may now be obtained from the poly- 
nomial fit 


r c = r ( V p c V 


and the equilibrium sound speed obtained from the relation 

r -j H 

i r c p c 1 
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For a Cartesian system, Ar in Equations (53), (54), (55) and 
(56) is replaced by Ax. The y momentum equation replaces 
Equation (57) and takes the form 


pu ( v x + H v z ) = - pvv y - p y 


(69) 


Noting that. 





(70) 


the expression for the'change in cross flow velocity along the 
streamline projection CD takes the form 


V C - V D ’ cos 


-X + tan if) tan ip q., 
pq \ y 


. !!*_ 
2 , 

COS 1 1 > 


Ax 


Properties at D are obtained using the relation 


(71) 


f 

D 



tan ij) 
icos <J> 



( 72 ) 


All other equations and remarks pertinent to the line source 
system also apply to the Cartesian system. 
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For the cylindrical system, the C point is located using 
Equation (73) 


r C = P D + tan ^ AX) 


( 73 ) 


Points A and B are located using the relation 


r A ■ r c - iXi h ix 

B B 

The pressure p q and flow deflection angle p ^ in the reference 
plane are determined by solving the linear relations along AC 
and BC with Ar replaced by Ax. 


The 0-momentum eouation may be written along the streamline 
projection in the reference plane, yielding 


v, 


P D r D q D C0S *0 


[ D . + pvv + pvq sin $1 Ax (75) 

6 W 


All streamline relations are valid along the actual streamline 
D*C, where properties at D* are obtained using the relation 


The streamline relations (Equations 62-68) described tor the 
line source system are also valid for this system. 


( 76 ) 
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E. Upper or Lower Boundary Point Calculation For a 
line source system*, the upper or lower boundaries are pre- 
scribed by equations of the form 

z = f(r,e) (77) 


where f is everywhere continuous and differentiable. Referring 
to Figure (16), which depicts an upper wall calculation, the 
appropriate boundary condition at point C is 


where 


and 


7 . n = 0 


[l j A S\ A 

v = u V + V i 6 + w i z 


Tf 


A 




A 



which yields 


(78) 


w 

u 

C 




where 

w 

— = udn w 

u Y C 

C 


(79) 


( 80 ) 


*C is the grid point IMAX(J), J at station r=r 2 and D is the grid 
point IMAX(J),J at station r=ri for an upper wall calculation. 
For a lower wall calculation the corresponding points would be 
the I,J points at the stations ^ and rj. 
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(J>j, being the flow deflection required for the characteristic 
calculation. The wall calculation requires the following 
iterative procedure: 

(1) A value for (^) is assumed, for examole the 

C 

value at point 0. 


(2) A value for <b r is obtained from Equation (80) 
and hence n can then be obtained using the 
compatibility relation along AC, (Equation 
55) with <J> C prescribed.* 


(3) The normal momentum equation (Equation 59) is 
applied along C-D yielding the cross velocity 

V C’ 


(4) The actual streamline massing through C is 

traced back locating the point D* on the ve- 


hicle wall and application of the pressure- 
density relation (Equation 62) (holding F 
constant) and the Bernoulli relation (Equation 


64) yields the total velocity v at C 


v 


2 





V 


* 

c 


*for a lower wall calculation the compati bi Vi Yy _ "relation along 
8C (Equation 56) would be employed. 


( 81 ) 
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(5) v and v* are checked for consistency and if 

they differ by more than a specified tolerance 

a new value of (-) is assumed and the process 

u C 

repeated until convergence is obtained. After 

two values of {— ) have been assumed, further 

u C 

values may be obtained by a standard linear 
error extrapolation procedure. 

For a Cartesian system, the upper or lower boundaries are pre- 
scribed by equations of the form 

z = f(x,y) (82) 

the velocity V is expressed by 

V = u i x + v i y + w i z (83) 


and the normal to the boundary is given by 


n 


f x i x + fy iy 


A 



(84) 


A 

Hence the boundary condition v.n = 0 yields 





( 85 ) 


The iterative procedure is the same as that described for the 
line source system. 
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For a cylindrical system, the upper or lower boundaries are 
prescribed by equations of the form 

r = f (9 »x) (86} 

The velocity V is expressed by 

V = u i x + v i e + w i r (37) 

and the normal to the boundary is given by 

A A f A A 

" = f x ix + ~ h - *r (88) 

Hence the boundary condition V n = 0 yields 

(it) = f x + (I) li (as) 

u c u c r 

The iterative procedure is again the same as that for the line 
source system. 

F. Sidewall Calculation For the internal flow 
analysis of Vehicle I, the sidewalls are calculated using a lo- 
cal cartesian reference plane system, the reference planes being 
the planes z = constant, as indicated in Figures (5) and (6). 
Additionally, for ail external flow calculations a cortion of 
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the vehicle undersurface can be interpreted as a sidewall in 
a local cartesian (y,z) system as indicated in Figure 
( 10 ). 

The sidewall calculation is geometrically complicated when a 
line source system is employed for the calculation of the other 
reference planes since local transf ormati ons must be made be- 
tween the cartesian and line source systems. This situation 
is depicted in Figure (17). For the internal flow sidewall 
calculation of Vehicle II, a local reference plane system is 
used where the reference planes are the surfaces r = constant, 
as indicated in Figure (7). The sidewall calculation is per- 
formed after all interior mesh points have already been calcu- 
lated. In performing this calculation it would be expeditous 
to use reference planes normal to the wall, but this is ruled 
out by the geometric complexities involved in obtaining cross 
derivatives for arbitrarily oriented planes. 

(1) Line Source or Cartesian System: Referring tc 

Figure (17), point C is located as the intersection of the 
sidewall with the plane z = constant at the station being calcu- 
lated. The sidewall is specified. by an equation of the form 


y - f (x,z) = constant 


(90) 
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The normal is then given by 


A /V 


n = [-f x i + j - f z k] /(I + f x 2 + f z 2 )^ 


(91) 


At point C, the wall boundary conditions is 


( v • n) r = 0 


(92) 


which may be written 


(-) 

u' 


f x + (W) f. 


(93) 


The streamline differential equations are 


dx _ _ dz_ 

u v w 


(94) 


Which in differential form become 


yc-yp 




(95) 


and 


*c- x d’ 


■* 1 ‘ft,* + 


( 96 ) 
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The z momentum equation, written along CD, takes the form 

as CD 

The characteristic compatibility relation (in the plane z=const) 
takes the form 

pu 2 d (^-) ± 8 dp * F* dr (98) 

where . 

F 1 = A 1 ( pw u z ) - pwv z + [A 1 - v ] (99) 

( “W p^ - a 2 pw£) 
a 

The properties at C are obtained by the following iterative 
procedure: 

Assume a value for (w/u) c ; for example, the value at the 
previous wall point E. Since f x and f z are known at C, Equa- 
tion (93) yields (v) for this guess. The location of the D* 

u c 

point may be obtained using Equations (95), (96) and the fact 
that this point has a radius r (or x is specified) and lies on 
the sidewall. It is obtained by interpolating between sidewall 
grid points. The D point (projection of D* onto plane z=constant) 
is then obtained by setting z Q = z c , y D = y D * and x D * x Q *. 


w c * w D 


SL + 
pq 


WW; 
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The A point is located using the characteristic equation on 
the plane z = z^ 


y C - *A 
X C ‘ X A 


2 

uv±a g 
U 2 - a 2 


A 


+ 



( 100 ) 


2 2 k 

and knowledge that r ft * r = (x ft + x A ) (or that x A is specified 
for a Cartesian system) and z^ =» z^, by interpolating between 
grid points on the plane z = z c . Obtaining properties at A, 
the compatibility relation (98) can be applied to yield p c for 
the assumed value (v/u)q. Streamline relations applied along 
D*C yield and v^.. The system can be checked for consistency 
as follows: 


V C’ (?> c and <5> c 1 u t 1 '( (101) 

* 

Test if w c = w^. to within a specified tolerance. If not, perturb 
the value of (w/u)^ assumed, repeat the above calculations and 

* 

once two values of (w/u) c and the corresponding errors (w c - w c ) 
have been obtained, convergence can be obtained by a linear 
extrapolation of the error to zero. 

Upon convergence, the remaining properties ( h c , . $ c » r c and a c ) 
can be obtained as described previously. 
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(2) Cylindrical System: The sidewall network for 

Vehicle II is depicted in Figure (18), This requires working 
in reference planes r = constant. The sidewall is specified 
by the equation 

e = f(r.z) 

The normal is then given by 

" = <- f z 1 z + ~ i e " f r + “2 + f r> 

A 

hence at point C, the boundary condition vn*0 yields 


<U>, = r c + ( U> ^ r 


The streamline differential equations are 


dr rde dx 
w” - v = u 


which in differential form become 


9 C ~ B P 
Ax 


(V) + (V) 

. U Q* U ! 


'K + r C ) 


r C - r D 
Ax 


1 

7 


(W\ + (E) 

u D* V C 


( 102 ) 


(103) 


(104) 


(105) 


( 106 ) 
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The characteristic directions, on the surface r = constant, 
are given by 

2 

qv cos 4>/a ±B 
(q/a) 2 cos 2 <j)-l 

where 

2 2 2 

q c = v c + u £ 

and 

= ( q/a ) 2 - 1 

Along these characteristics, the compa ti bi 1 i ty relations take 
the form 


de = 1 
dx r 


d (^.) ± 3 dp = F* dx 

where F 1 is given by 

F* = p ( - w v r + X ± w u r ) - (X ± u - v) 

( pWy* + PfO 


The normal or w-monentum equation may be written along a stream- 
line projection in the reference surface r = constant. 


W C = W D 



2L + w w r 
P r 



(107) 


(108) 


(109) 


AS 


( 110 ) 
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where 

h 

AS = (r 2 (e c - 9 0 ) + Ax 2 ) 

The logic for the sidewall calculations proceeds as follows: 
assume a value for (u/u)q. Using Equation (104), yields 
(v/u)^. The streamline point D* may be located by Equations 
(105) and (106) in an iterative manner. Using Equation (107) 
to locate the A point, we may apply Equation (108) along AC 
yielding p^, and then apply Equation (110) to determine w^. 
Using streamline relations along D*C yields Pq and v^. The 
system is then checked for consistency; since vq, (w/u)q and 

1c 

( v/u) c yield an independent value w c which is compared with 
the value of Wq obtained from the normal momentum equation. 
The same extrapolation method can be used for convergence 
and the remaining properties (h c , 3> c , r c and a c ) may be ob- 
tai ned . 


G . Internal Corner Flow Calculation 


(1) Line Source System: The corners are de- 

fined by the intersection of the upper or lower nozzle walls 
with the sidewalls. The upper or lower wall is specified by 
an equation of the form 


z - f(r,e) 


(111) 
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while the sidewall is specified by 
y = g(x,z) 

The transformation from a Cartesian to a line source system is 
simply 


x = r cos 0 
y = r s i n 8 

z ~ z 

With u, v and w being the velocity components in the line source 
system (in the r, 6 and z directions resoecti vely ) and u\ v and w 
being the components in the Cartesian system the velocity trans- 
formation is 

IT = u cos 9 - v sin 6 
v = u sin 8 + v cos 0 
w = w 


At the corner (point C of Figure 19) the boundary condition 

_ A 

v . n = .0 applied to the upper (or lower) surface yields 





Applying this to the sidewalls yields 



flx + '<=> 9 Z 


( 112 ) 


(113) 


(114) 


(115) 


(116) 
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The vertical reference plane 0 = 0 C nay be intersected with the 
initial plane as shown in Figure (18). 


Applying the velocity transformation (Equation 114) to Equation 
(116) we obtain 


d\ g x + (tt) 9 z tan 9 
u (1 + g x tan 9) 


(117) 


Nov;, Equations (116) and (117) can be solved for (jj) yielding 


{-) = 
u 


f r (1 + g x tan 9 ) + f 9 (g x - tan 9) 
(1 + g x tan 8 - g z p) 


(118) 


Then 4 > q = tan * (]j) and the compatibility relation (Equation 55 
or 56) can be applied along the projection of the bi -character! s ti c 
onto the vertical reference plane ( A 2 C ) yielding pq^. The point 
A 2 is located by a double interpolation. 


The horizontal reference plane z = z c may be intersected with the 
initial plane as shown in Figure (18). Transforming Equation 
( 1 1 5 X to the Cartesian system, we obtain 

(=) = cos 9 f = sin 0 £1 + (-) [sin 9 f + cos 0 £§.} (119) 

u r r xr r r 

Solving Equation (116) and (119) for (~) we obtain 

u 
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v 

u 


g + [cos ef -sin0f6]q 
*x r y— 2 3 z 


1 - [sin e f r + cos e — ] g • 


Applying the compatibility relation (Equation 98) along the 
projection of the bi-characteristic onto the horizontal refer- 
ence plane (AjC), we obtain an alternate value of the pressure 
Pq^, where the point Aj is again obtained by a double interpola 
tion. The pressure p^ is then obtained as a weighted average 
of the two pressures p^ and p£^» wei 9^ 1n 9 function being 
determined by the relative wave strengths along the two char- 
acteristic projections AjC and A 2 C. 


Since a sharp corner is a streamline of the flow field, the 
streamline relations may be directly applied along CD yielding 
p c , v c , h c and Knowing the velocity v c and the flow angles, 

we obtain the velocity components u, v and w or u, v and w. 

2. Cartesian System: The upper or lower wall is 
specified by an equation of the form 

z = f(x,y) 

No velocity transformations are required in this calculation. 

The boundary condition v • n applied to the upper (or lower) 
boundary at C yields 


( 120 ) 


( 121 ) 


( 122 ) 
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while applied to the sidewall yields 


(-) =o 01) a 
K u ' - x v u - z 


Solving for (^-) one obtains 


( 123 ) 


(-) = 
1 u' 


9x ^y 
' 9z f v 


(124) 


The vertical reference plane y=y c is intersected with the initial 
plane and the characteristic compatibility relation (Equation 55 
or 56) is applied along the projection of the bi -characteristic 
yielding p^* since <f > q = tan * (tj) is known. 


Similarly, the horizontal reference plane z=z r is intersected 
with the initial olane and the characteristic compatibility re- 
lation (Equation 98) is aoplied along it yielding an alternate 
value of p^ where (-■) is given by 


(1) g x + f x Sz 

U " 1 " f y.9z 


(125) 


The pressure is again a weighted average of p r and based on 

U 1 ^2 

the relative wave strengths along the two characteristic pro- 
j ecti ons . 


3. Cylindrical System: The upoer or lower wall is 

specified by an eouation of the form 
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r = f (x,6) 

while the sidewall is specified by 
6 ■ g(x,r) 


At the corner, the boundary condition v • n apolied to the upper 
(or lower) surface yields 


(?) - f x ♦ <*> f£- 

and applied to the sidewalls yields 

(J) = 9x + (ff> 9r 


The reference plane 9=9^ passing through the corner point is 
intersected with the initial plane and the characteristic com- 
patibility equation (Equation 55 or 56) is applied along the 
projection of the bi -charachteri sti c onto this reference plane. 
This yields o since <j> is known, being obtained from the re- 

'4 4 

1 a t i o n 



tan" 1 (]j) - tan" 1 


f -\ 

X x r 

l - g r 

r / 


( 126 ) 


(127) 


(128) 


(129) 


(1.30) 


Similarly the surface r=r^ passing through the corner is inter- 
sected with the initial plane and an alternate value of the 
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corner pressure p c is obtained using the compatibility rela- 
tion (Equation 108) along the projection of the bi -character!* sti c 
onto this reference surface, where (v/u) is given by 



* f x Sr 
i f 9 „ 


(131) 


The pressure is then again determined by a weighted average. 

H . Flow Field Discontinuities 

The numerical program developed can perform three dimensional 
calculations for discontinuity surfaces propagating into a 
nonuniform stream. The use of the reference plane technique 
can greatly simplify the logic by using a two dimensional type 
of procedure in the reference planes. While this procedure re- 
sembles its strictly two dimensional counterpart it should be 
emphasized that the full three dimensional Hugoniot relations 
are satisfied at a shock surface, and full three dimensional 
boundary conditions are applied at a contact surface. 

Consistent with the problem considered, only discontinuities 
which propagate essentially in references planes can be analyzed. 
This implies that discontinuities formed from sidewalls or large 
cross-wise flow gradients cannot be analyzed. Equilibrium chemis- 
try is assumed to prevail upstream and downstream of a disconti- 
nuity.. 
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(1) Shock Surface: A typical local shock sur- 

face orientation is shown in Figure (20). The local shock sur- 
face is specified by an ortho-normal triad of vectors consisting 
of the normal to the surface and two tangent vectors to the sur- 
face. With respect to the reference axes, a tangent direction 
t with cosine director 8 is defined by the shock and reference 

A 

plane intersection. The tangent £ with director (a) is normal 

A A 

to t and the normal to the shock n is n = t x £, or in terms of 
a and 3 


n = - cosasing i r - sina i 0 + cosacosg i z 
The velocity vector in terms of the shock oriented coordinates 

'A ~ A "'A 

is v = u n + v t t + v A l. The local Hugoniot relations are: 


Pi Uj 

“ .92 u 2 

continuity 

~2 

~2 


Pi + Pl u 2 

= P 2 + P 2 u 2 

normal momentum 

Vt l = 

Vt 2 

/V 

t momentum 

V *1 " 

v £ 2 

A 

£ momentum 

H = constant = h+% V 2 

energy 


(132) 

(133) 

(134) 

(135) 

(136) 

(137) 


p = p (P,h,$) 


state 


(138) 
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Combining and rewriting (133) and (134) in terms of the refer- 
ence plane parameters we obtain, 


pj u 1 (u^-u 2 ) = p 2 " P i continuity & n (139) 

momentum 

( u = v s i na+qcosas i n( 8-<f> ) ) 

q-^ cos(6-<j)j) = q£ coS'( 6-<J> 2 ) . t momentum (140) 

/V 

qj (-sinasin(e-<|> 1 ) + cosa tan ij)j) l momentum (141) 

= d2 ( -s i n aS i n{ 8-02) + cosa tan ^ 2 ) 

2 

H2 = h 2 + H (q 2 /cos tp 2 ) = constant energy (142) 

P2 = p( p 2» h 2’ state (143) 


A typical shock wave calculation would be performed as shown in 
Figure (21). A new shock location (point C) is determined in 
each reference plane using an assumed value of the cosine direc- 
tor 6 C * A new value of the cosine director a c can now be calcu- 
lated from geometric considerations using the newly calculated 
shock locations, and the relation tana = tana' cos8. Since prop- 
erties are nonuniforrr upstream of the shock wave, a reference 
plane characteristic calculation is performed on the upstream side 
of the wave (AiC,BjC) and local flow properties determined at Cj. 
Using this data and the a and 8 cosine directors, the Hugoniot 
relations (Equations 139 through 143) are solved in each reference 
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plane for the downstream flow conditions. Satisfying the 
Hugoniot relations involves an iteration procedure since chemi- 
cal equilibrium is assumed upstream and downstream of the shock 
wave. This particular procedure is analogous to that employed 
in two dimensional flow calculations. Using the calculated 
Hugoniot flow deflection in each reference plane a value of pres- 
sure is calculated on the downstream side of the wave ( C 2 ) using 
the compatibility relation along ( C 2 - A 2 ) - This pressure is com- 
pared to the value calculated from the jump relations. If the 
difference in pressures exceeds a specified tolerance, a new 
shock wave location is computed using a new value of the cosine 
director & c for the reference plane being calculated until con- 
vergence is achieved. This process is repeated in all reference 
planes. After convergence is obtained in all planes, improved 
crosswise directors a c are computed using the converged shock lo- 
cations and the procedure may be repeated until two subsequent 
values of a in a given reference plane agree to within a speci- 
fied tolerance. 


(2) Contact Surface: The numerical grid for the 

computation of a three dimensional contact surface is depicted in 
Figure (22). Let (1) and (2) denote conditions below and above 
the contact surface respectively. Properties are to be determined 
at the point C in the reference plane 9 = constant. The contact 
surface cuts this reference plane along the line C-D, the local 
angle of the cut making an angle £ with respect to the r direction. 
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The local geometry of the contact surface is identical to that 
of the shock surface as depicted in Figure (20). 


There are two streamlines passing through the point C: C^Dj* on 
the lower surface and C 2 D 2 * on the upper surface. These stream 
lines project onto the reference plane as C 1 D ^ and C 2 D 2 re- 
specti vely . 


The appropriate boundary conditions at point C are: 




where 


A A A A 

n = - cos a sin 6 i r - sin a i 0 + cos a cos 6 i z 


tan a = tan a' cos 3 


and 




(144) 

(145) 

(146) 

(147) 


( 143 ) 
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The numerical solution proceeds as follows: 

(a) In each reference plane e=constant assume a 
value for B c » for example B c = Bg* The point 
C can then be located and values of a c can be 
computed using the geometric location of the C 
points in all reference planes to compute a', 
and Equation (147) to determine a. 


(b) A local iteration loop is required on either 
side of the contact to satisfy the following 
system of equations for the assumed value of 8 
Equations (144) or (145) take the form 


sin (8-4>„ ) + tan a r 
C 1 C 
7 


tan ^ 
7 


0 


(149) 


The compatibility relations along AC or BC take 
the form 



The Bernoulli, relation (Equation 36) applied 
along CjDj* and C 2 D 2 * yields V c ^ and V^. 
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The normal momentum equation applied along 
CjDj and C 2 D 2 yields v c ^ and v C2> 

The cross flow angle i|> is obtained from the 
relation = sin - * (v/V) (151) 

(c) The iteration proceeds as follows: A value of 

Beta ( B c 2 ) at C 2 is assumed and a pressure 
(P C1 - P C2 ) is guessed. For this pressure and 
Beta, Equations (150), (151), and (36) yield all 
properties on either side of the contact. The 
boundary condition (149) on side 2 then is itera- 
ted until a consistent set of values B C2> p c , 

<$> r , ^) r and a are obtained. 
c 2 c 2 

The boundary condition (149) on side 1 then yields 
a value for B c ^ If this does not agree with B C2 
to within a specified tolerance a new value of 
pressure is guessed. From a linear error test a 
new value of Beta is obtained and the process re- 
peated until Equations (144) through (149) are 
satisfi ed . 



TR 166 - V o 1 . I Page 77 

I . External Flow Interaction 

As a result of the cowl length being considerably foreshortened 
a substantial amount of nozzle under-expansion can be expected. 
The resulting wave patterns can at certain flight conditions in 
fluence the pressure distribution on the vehicle undersurface 
and hence the thrust, lift and pitching moment. Figure (23) 
illustrates the nature of the waves produced in a typical refer 
ence plane. In some cases, as with end module configurations, 
significant sidewise under-expansion may also occur. The re- 
sulting waves are no longer quasi two dimensional but fully 
three dimensional in nature. 

(1) Cowl Under-expansion Conditions At the 
cowl edge the flow phenomena in planes normal to the edge is 
locally two dimensional. Thus, in these planes the cowl under 
expansion shock and Prandtl -Meyer expansion can be computed 
from two dimensional considerations. It is assumed in these 
calculations that the chemistry is in equilibrium. Figure (24) 
depicts this locally oriented system. For a given cowl shape, 
data from the reference system may be transformed to the local 
coordinate plane. The calculations are performed in the local 
system and then transformed back to the original reference sys- 


tem . 



TYPICAL REFERENCE PLANE 



FIGURE 23, UNDER-EXPANSION INTERACTION 
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In Figure (25) let Ji be defined as the local unit tangent vector 
to the cowl edge. Then for a Cartesian (x*y,z) or cylindrical 
reference system (x,0,r): 


cosa 

A 

V.e 

+ sina i z>r 



sine 

A 

\ - 

cose 

Sina ? y,e 

+ cosB cosa 

A 

Vr 

cose 

A 

- 

sine 

A 

sina i., a 

y » w 

+ sine cosa 

A 

1 z ,r 


where n is the unit normal to the wave and t is the second tan- 
gent to the wave surface. Now: 

oil = Arctan or I lL) 

x dy r de 

6| = Arctan (^ or ~£) 

y,e dx dx 

and 

tang * cosa tang 

The total velocity is defined as 

s q(cos<|>i +tani|)i^ Q + si n4>i ^ ^) = u n + v A + w t< 

Then : 


(152a) 

(152b) 

(152c) 

(153a) 

(153b) 

(154) 

(155) 
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A 



FIGURE 25. EXTERNAL FLOW INTERACTION-LOCAL AND REFERENCE 
PLANE ORIENTATION 
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u = - q ( - s i n 6 c o s 4> - cosj3 sina t a n4> + cos6 cosa sin$) (156a) 

v = q(tamjj cosa + s l* n4> sina) (156b) 

w = q(cos£ cos<p - sinB sina' tamp + sinfi sin<}> cosa) (156c) 

For locally two dimensional flow \f = constant. Then for the ex- 
pansion region the equations to be solved are: 


d£ 

+ i q 2 = 0 

(157a) 

P 

2 


P 

II 

( i 57b ) 

dfr 

a. 

= - u 

(157c) 

de 

.2 


h + 

= constant 

( 1 5 7 d ) 


2 

^ 2 2 0,2 2 ^ 2 

where: u = ( q -w ) = ( r £) , q = V and T is the equilib- 

P 

rium isentropic exponent. 

For small steps Ap in pressure the solution for r constant over 
the integration step 

i 2 = - yf±T (k 2 -k 1 ) + i l 


(158) 
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where 

k * Arc tan j — 

\/i- (w/c )2 

and 

c 2 . iii Sf 2 ♦ a 2 

r-i i 

Thus the expansion region may be solved in a step by step manner 
using the above solution. The numerical iteration scheme is as 
f ol 1 ows . 

Assume a shock strength Beta($). Using the local Hugoniot solu- 
tion previously described compute all flow properties downstream 
of the shock using the cowl external data. From the pressure 
downstream of the shock and the cowl internal data compute the 
Prandtl -Meyer turn using the solution described above. If the 
deflection angles from the shock and Prandtl -Meyer turns do not 
agree to within a specified tolerance a new shock strength is 
assumed and the process repeated. 

After the above solution is obtained the local properties are 
transformed back to the original reference systems. That is: 



qcos<f> = u = u sing + ft cosg 


(159a) 
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q t a nip = v = u cose sina r v cosa 

- 5? sine sina 

(159b) 

qs i n$ = w * u cose cosa + v sina 

+ w sine cosa 

(159c) 

2 2 2 
q = u +w 


( 159d ) 

tans = tane/cosa 


( 159e ) 

(2) External Corner-End Module 

As long as the 

cross 


flow in thevicinity of the external corner region remains con- 
tinuous a local r , e system can be used to represent the angular 
property variation at a given marching step Figures (26a) and 
(26b) depict the geometry of the reference system for an end 
module configuration. 


For any property the crosswise solution may be represented as 
a polynomial 


f = a + bs +--+s n 

and 

r = cn (e ) + g(e) 


(160) 


where "$" is a distance along a constant n surface. Since only 
the existence of the first derivative is required f and r may 
be approximated by 


f = a + bs 


and 

r = c(n Q + nje) + g(e) 

In region I we define (see Figure 26c) 


( 161 ) 
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n I * ^contact - 1 2 sav> (l62a> 

g(e) « 0 

For region II we define 

n II * ^contact - z shock ) (162b) 

*’ r contact 

and for region III we define 

n III ■ (z max ' z shock> <16Zc) 

9 (e) - ""shock 

The above technique allows the external corner region to be 
scaled in a manner consistent with the growth of the shock and 
contact surfaces, and analytically represented at each marching 
step . 


J . Swept Throat 

Due to vehicle aerodynamic and combustor design considerations 
it may be advantageous to introduce sweepback into the flow 
field upstream of the nozzle throat. The program described in 
this report is capable of accepting initial data along a 
swept initial surface. The following limitations apply to the 
sweepback option. 

(1) The flow normal to the sweepback is super- 

sonic. 
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(2) The sweepback is linear; that is, the cross- 
wise location (y or 9) of the sweepback 
line on the upper or lower surfaces is a 
linear function of the marching distance. 

(3) The sweepback lines lie on cylindrical or 
Cartesian surfaces. 

(4) The sweep of the upper and lower surfaces 
i s i denti cal . 

Figure (27) illustrates the geometry associated with a typical 
swept throat as defined in this program. 

The numerical procedure is as follows: all data along the 

initial swept surface is stored in locations J = 1, 3, 4 , 

Initially calculations are carried out only in the J=1 plane 
with the necessary data for derivatives being interpolated be- 
tween the first and second (stored as third) plane on the ini- 
tial surface. This procedure is carried out until the marching 
step has reached the second reference plane. The marching con- 
tinues with both the first and second planes being calculated, 
the data necessary for calculating derivatives on the second 
plane are interpolated between the second and third (stored as 
fourth) planes. Thus as the marching proceeds each subsequent 
reference plane is added to the calculation at the proper loca- 
tion. This continues until the swept region has been completed, 
whereupon the calculation automatically reverts to a normal 
marching procedure. 
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FIGURE 27. SWEPT THROAT 
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K. Discontinuities at Sidewalls 

Implicit in the use of the reference plane technique is the 

assumption of a continuous crosswise flowjthus discontinuities 

formed at sidewalls or by large cross flow gradients cannot be 

analyzed. Only discontinuities which primarily propagate in 

the reference planes are considered. Consistent with this 

assumption no reflections of discontinuities at sidewalls are 

permitted. This is expressed geometrically as n * n = 0. 

s w 

Where n $ is the normal to the discontinuity and n w is the 
normal to the sidewall. Figure (28) illustrates the geometry 
of the sidewall and discontinuity intersection. 


Assume that the discontinuity may be described by an equation 
of the form 


z = ^ ( X , y ) 


( 163a) 


n 

s 


i 

z 



/ \ 



and the sidewall by an equation of the form 

y = g ( y » z ) 



(163b) 


(164a) 


(164b) 




FIGURE 28. DISCONTINUITY AT SIDEWALL 
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Then n *n =0 yields 
$ w 


f x 9x - f y - 9 Z - o 


From Figure (28) along (1-2) we may write 


Z 2 ~ 1 


1 + f -~ 1 • <V‘l> * fj ^L . <vV 


and along (3-2) we may write 


Z 2 = Z 3 


+ '*£21 . ( V , 3 , 


Equations ( 164a ), { 165 ) , ( 166) and (167) are four equations in 

four unknowns y , z ' t f , f . 

X 2 y 2 

The numerical solution proceeds as follows: 


(1) Assume z^. 


(2) Using Equation (164a) compute y„, g and g 

^ Xa i 


(3) From Equation (167) compute f and from Equa^ 

y 2 

tion (166) compute f x . 


(4) Equation (165) is checked for consistency. 


( 165 ) 


(166) 


(167) 
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If Equation (165) is not satisfied to within a given tolerance 
then a new value of z 2 is guessed and the process repeated 
until convergence is achieved. 


L . Thrust, Lift, Pitching Moment 


The following definitions are used in this report for thrust, 
lift, pitching moment 


T 



V dA n 


(168) 


L 



V dA n 


(169) 



(170) 


Figure (29) gives the orientation of the vectors with respect 
to the vehicle. Internally the integrals range over all the 
vehicle surface areas. Externally they range over the complete 
vehicle undersurface as defined by the bounding stream surface. 
For external flow with wraparound, val ues of these parameters 
are computed upto the contact surface and shock surface, in 
addition to the values at the bounding stream surface. 



V- ENGINE 
MODULE 



FIGURE 29. THRUST, LIFT, MOMENT 
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Based on throat data the program computes at the throat 
station values for the stream thrust in both the z and x 
directions. These quantities are defined as 



T s* • (pV-S) dA n p-pj 1 x .dA n 

A- 



sz 



V • i (pV* n) dA + p-p ) i »dA 

c n 00 z n 



M. Integral Correction Factors 

Numerical solutions of the differential equations of motion 
generate errors associated with machine round-off procedures 
and truncation errors corresponding to the numerical algorithm 
employed. These errors result in global conservation of mass 
and energy not being satisfied exactly. The subject computer 
program can as a user option limit the propagation of these 
errors by providing at arbitrary calculation steps integral 
corrections for mass flow and energy. These correction pro- 
cedures are described below. 

Conservation of mass and total energy require that 
M = 



n) dA n = constant 


(171) 

(172) 


(173) 
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» JJ * (h + \ V 2 ) (pV • n ) d^ n = 


constant 


Referring to Figure (30) the integrals may be represented by 
the summations 

M = l p u AA 


i .J 


and 


H = i (h + I V 2 ) ( PU ) aA 

i,j 2 X 

where the barred quantities refer to averages over a numerical 
mesh. The average of any property f over a mesh interval can 
be defined as 


Z . ./f ds 

f 


where s is the line distance surrounding an elemental area dA. 


The integral correction factors for mass and energy can be de- 
fined by noting that 


n * M + e l 



C 1 / / < h + 7 V > ^ V * n) dA n * H + e 2 


(174) 


(175) 


(176) 


( 177 ) 
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Z 


- Y 



FIGURE 30. ELEMENTAL AREA ASSOCIATED WITH GRID SIZE 
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where ej and c 2 are errors In mass and energy. 

Then 


, . M C - • M 

M+e ^ ’ 2 M+c 2 


(178) 

With and C 2 defined the following system of 

solved for the corrected properties 

equations are 


* * T 

p /p * constant 


(179a) 

p*/p*= RT* 


(179b) 

h* = h(p*, <r>, T ) 


(179c) 


p* V* = Cj P V 

(h + \ V 2 ) = C 2 (h + \ V 2 ) 

Thus, properties of each point in the flow field are corrected 
uniformly for the integral correction factors for mass (C^ 
and energy (C 2 ) . 


( 1 7 9 d ) 
( 179e ) 
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N. Second Order Procedures 


The numerical methods and analytic techniques presented in the previous 
sections were for the purposes of clarity and simplicity assumed to be 
first order in mesh size(Ax, Ay). However, the program is capable of 
second order computational accuracy with respect to both local and 
global calculations. Local refers to the computation of a single mesh 
point in a given reference plane using a fixed value of the cross deriva- 
tive. Global refers to the computation of all the mesh points in the 
flow field with updated values for the cross derivatives. 

1. Interior Point Calculations 

Local Iterations - The basis for establishing a local 
second order numerical characteristic type calculation is well known (see 
Reference 2) and will briefly be described below. 

Referring to Figure (12a) a typical interior point calculation for the line 
source system is performed as follows. The point C is located by projecting 
the quasi -streamline a distance Ar downstream. That is, in difference form 
Equation (22) is 

z c = z D + (o tan $ D + $ tan <j> c ) Ar 

where a and $ are the local averaging coefficients with values (a « 1, $ * 0 
or (a = h» 6 s h). The initial computation (first order) is made with 

a - 1, & = 0. Points A and B are located by shooting back, characteristic 


( 180 ) 
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projections from point C, and locating these between points 1-1,1 and 
1+1,1 respectively, by an iteration procedure. Using Equation (21) in 
difference form: 

Z A,B = Z C “ X A + ^ Ar ’ 

B 

All required properties can be obtained at A and B including cross deriva- 
tives by a linear interpolation procedure. 

The pressure p^ and flow deflection angle in the reference plane can 
be determined by solving the compatibility relations along AC and BC. 

(5 (pq\ + B (pq 2 ) c ) <* c - + A > MS B a + h c ) (P C - P A ) 

(S (pq 2 ) B + H (pq 2 ) c > (* c ' *b' ' 6 b + ® B c ) tp c ‘ p e) 


* f b Ar 


where for a purely local iteration F* are evaluated on the initial value 
surface. The cross velocity v^, may §e obtained from Equation (59): 


[ / a \ ■ / % \ p( s ) 

V r. ~ V n • i ( y* rn c A* ' P 


X D [ cos $ 


T COS ^ 


where F^ is evaluated only on the initial plane for a purely local itera- 


(181) 


(182a) 


(182b) 


( 183 ) 
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tion. 


Properties at 


may be obtained from the relation 





•) A0 


which from geometry becomes 



tan 0 x . / f tan q , * Ar 

! v 0 r cos <t> ' D v e r cos <j>'D 2 


(184a) 


(184b) 


The pressure density relation (13) may be integrated along 0*C holding r 
constant for a small integration step, yielding as in Equation (62) 

P C * p 6* 

The Bernoulli relation, Eqution (64), then yields the total velocity v c< 
Conservation of total enthalpy, Equation (65), yields h c and for in- 
viscid flows from Equation (66) 



* c * V 

The isentropic exponent may be obtained from a polynomial curve fit 
r c = r (h c> p c> $ c ) 

and the sound speed from Equation (68) 
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For a purely first order calculation (a ■ 1, £ = 0) the computation would 

now be complete. However, for a local second order calculation, the above 

process would be repeated with (a * 6 = h) until the newly calculated 

pressure p^, converged to the previously calculated pressure p^ v “^ to 

+ ( s ) 

within a given tolerance. Note that in a purely local iteration F“, F v 
are held constant at their Initial values. 


Global Iteration - After completing all calculations if 
overall second order accuracy is required the above process is repeated 
with updated values of the cross derivative functions. That is we may 
write: 


F AC * a ' F A + b ‘ F C 
BC B 


(185) 




+ 



(186) 


where a and b are the global averaging coefficients with values, a = 1, 
b = 0, for a first order global and, a ® h> b = for a second order 
global calculation. The cross derivative functions F£, F^ are evaluated 
on the new data surface (^ + Ar) using the converged properties previously 
computed from local iterations. These functions are held constant while 
new local iterations are performed with a = h» 6 = This process requires 
storing the cross derivatives at the new data surface in addition to the 
flow field properties. 



TR 166 - Vol . I 


Page 104 


The program thus provides the user with a choice of computational ac- 
curacy which can be tailored to the complexity of the flow field. That 
is, in some instances it may not be necessary to go beyond the local 
iteration procedure with consequent savings in computational time. Thus, 
computational accuracy from first order to fully second order may be chosen 
at the users discretion. 

2. Boundary Point Calculations 

All wall boundary calculations are formally identically to those described 
in the earlier sections of this report. Second order procedures are ap- 
plied to characteristic coefficients and cross derivative functions, as 
described above for an interior point calculation. The local iterate is 
the pressure p^. 

A fully second order (global ) shock and contact discontinuity calculation 
requires updating of both the cross derivative functions and the discon- 
tinuity surface cross anglejupon completion of all the local calculations. 
All local discontinuity computations remain formally identical to those 
described earlier in this report. However, the characteristic relations 
are modified as described above for locally a second order interior point 
calculation. For shock wave and contact calculations the local iterate 
is the surface discontinuity angle, Beta, in the reference plane. That 
is locally converged properties are obtained when Beta and Beta^ v_1 ^ 
agree to within a specified tolerance. 
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0 . Embedded Shock Waves 

The presence of compression waves in a supersonic flow field 
may, depending on the numerical scheme employed, lead to com- 
putational problems if these waves are not treated as finite 
discontinuities. The Hartree type grid employed in the numerical 
scheme does not follow characteristics. This property of the 
scheme renders it relatively insensitive to numerical diffi- 
culties if moderately weak compressions are encountered. 

In each reference plane a two dimensional type of procedure is 
employed to detect characteristic crossings. In each plane (J) 
using the initial data at points (I) and (1+1) the following 
equations are solved for a downstream crossing (if any) of waves 
of the same family (see Figure 31a). 

Z C ■ z a + ( V x a> 

z c * z b + tanU±M) b U c -x d ) 

The minimum downstream location of these crossings is determined 
and stored for later use. This procedure is repeated in each 
reference plane and an overall minimum crossing (i.e.« strongest 
wave) is determined. The initial data associated with this 
crossing will determine whether a strong compression wave is pre- 
sent at this point. Figure (31b) indicates the isentropic two 
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FIGURE 31a. CROSSING OF CHARACTERISTICS 
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dimensional flow deflections associated with waves, A($±y), 
propagating at various Mach numbers. 

A test is now made of the strength of the strongest wave de- 
tected with respect to a specified tolerance. If the wave 
strength, A($±y) exceeds this tolerance a shock wave strength 
is assumed to be of the average UtyJj and U±y)j + 1 . That is, 
s shock + (4* ± u) I + 1 )/2 between points (I) and (1+1). 

It should be noted that both grid spacing and wave strength 
tolerance determine the location of the shock formation. Once 
detected the program does not carry the shock wave as part of the 
internal calculation but does specify its location. 

P . Simulator Fluids 


The computer program described in this report contains provi- 
sions for the use of two gases other than equilibrium mixtures 
of Hydrogen-Ai r . As discussed in Volume IV of this report CH^ 
(Methane) and (Ethylene) can be used in cold flow exhaust 

nozzle tests to accurately simulate various equilibrium Hydrogen- 
Air thermodynamic and nozzle flow parameters. Curve fits of the 
thermodynamic properties of CH^ and have been taken from 

data given in Reference (2) of Volume IV. 

The computer program will accept CH^ or as the nozzle simula- 

tion fluid by simply specifying hypothetical fuel-air equivalence 
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ratios 4 = - 1 (CH^) or t » - 2 ^ C 2 H 4 ) for the ^ riterrial flow. 
For the external medium 4 » 0 specifies air 1 nequi 1 1 bri urn. 

All subsequent logical and numerical operations of the program 
remain unchanged. 
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IV. SAMPLE CALCULATIONS 

In order to demonstrate the program's capabilities, three 
sample calculations have been performed. 

For Case I, the flow field in a Vehicle I central module 
was computed. For the internal portion of the flow field, 
a line source reference plane system was employed as shown 
in Figures (32) and (33). Five reference planes (J=l, --5) 
were employed in this calculation with eleven data points on 
each reference plane. The radius of the initial station was 
equal to 27.5127. The external flow was taken as representa- 
tive of that on the surface of a 10 degree cone at an altitude 
of 165,000 ft. and Mach number of M=10. 

The streamline pattern and shock and contact shape are shown 
in Figure (34) for reference plane J = 1 , while axial surface 
pressure distributions in this reference plane are plotted in 
Figure (35). 

An analysis of the pressure distribution on the vehicle under- 
surface in Figure (35) indicates that the first expansion wave 
from the initial turning, reaches the undersurface somewhat 
downstream of the cowl exit. In addition, some reflected 
waves from the internal portion of the cowl appear to be cap- 
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FIGURE 34. CASE I OVERALL SHOCK, CONTACT AND STREAMLINE PATTERN - PLANE 1 £ 

n> 


co 


TR 166 



TR 166 - Vol. I 


Page H4 



X 


FIGURE 35, CASE I AXIAL SURFACE PRESSURE DISTRIBUTION 
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tured by the vehicle undersurface. The last waves appear to 
reach the undersurface at x=46, just short of the vehicle end 
It should be noted that no waves from the external flow inter 
action reach the vehicle undersurface. In Figure (36)* the 
Mach number and pressure profiles are depicted at the cowl 
end. 

The flow variables at the initial station (KOUNT-O), a sta- 
tion just downstream of the cowl end (K0UNT=55) and at the 
vehicle end (KOUNT-114) are tabulated in Tables la, lb and 
Ic respectively. The reference planes are denoted by J=l, 2 
etc., and the grid points in each reference plane by the index 
I: the reference planes in the internal flow are in a line 

source system and hence each reference plane has a constant 
6, while for the external portion of the flow a Cartesian 
system is employed and each reference plane has a constant 
value of Y as indicated. 

The value of the Z coordinate in each reference plane is tabu 
lated under Z, the pressure under P, the velocity component 
in the reference plane under Q, the flow angle in the re- 
ference plane under PHE, the cross flow angle i p under SI, the 
Mach number under M, the static enthalpy under H, the fuel to 
air equivalence ratio $ under PHI, the density p under RHO, 
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FIGURE 36. CASE I PRESSURE AND MACH NUMBER PROFILES AT COWL STATION 
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the equilibrium isentropic exponent r under GAM and the tem- 
perature under T. In the external flow data, the shock and 
contact are identified as points where two successive I 
Indices have the same value of Z: the shock has a jump in 

pressure for these indices while the contact has the same 
value of pressure. 

Case II demonstrates a fully second order calculation of an end 
module flow field. A simple two dimensional expansion of the in- 
ternal cowl surface was chosen, with a uniform external flow. This 
generated mild lateral pressure gradients at the cowl edge but 
significant corner gradients externally. 

The initial pressure was assumed to be 500 lb/ft at a Mach 
number of 2.21 and an equivalence ratio of 0.8. The external 
flow conditions were those of Case I. The lower surface and 
sidewalls are flat surfaces, while the upper wall has an ini- 
tial small radigs turn of 10° followed by a flat 10° surface 
to the cowl edge. The flow variables at the initial station 
(K0UNT=0), a station at the cowl end (K0UNT=70) and the last 
calculated station (K0UNT P I20) are tabulated in Tables 1 1 a , 

1 1 b and lie. For the external wraparound flow the tabulated 
values of Y on the rotated planes are measured from the last 
internal reference plane (J-JINT). That is Y = Z. 
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As shown in Figures (37) and ( 33 ), significant pressure gra- 
dients are induced in the external corner interaction region. 
Contract surface pressures equilibrate more quickly than those 
downstream of the shock surface. This is the result of the con- 
tract surface being a stream surface of the flow field. 

Case III demonstrates the programs' sweepback capability. The 
initial conditions are similar to those of Case I but are given 
along a swept initial surface. A cylindrical type of configura- 
tion was chosen having geometry similar to Case I with the initial 
surface swept back approximately 30 degrees. The flow variables 
at the initial station (K0UNT=1) and at a station downstream in 
the recompression region (K0UNT=110) just downstream of the end 
of the sweepback region are given in Tables Ilia and Illb. All 

surfaces are identical along lines (X-X ) equal to a constant. 

swp 

The sidewalls are assumed flat along a surface theta ( 0 ) equal to 
a constant, thus the expansions should be two dimensional along 
lines (X-X swp ). Figure (39) indicates that the program is calcu- 
lating correctly. Variations in pressure between reference planes 
that occur downstream of the expansion region are due to reflec- 
tions of waves off the internal symmetry plane. The flow in this 
region has become fully three dimensional. 
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FIGURE 38. CASE II SHOCK SHAPES AND PRESSURE DISTRI- 
BUTIONS AT VARIOUS AXI1L LOCATIONS 
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2,2l\p+00 6, 450^*07 fl OOOp.Ol 

2,2lfp,00 6,4bO E *n 7 0 t OOO F .Ol 

2,21\p*00 6,450 E ,i)7 e.OOOp-01 

2,21 i F , 00 8,450£+O7 e OOOp-Ol 

2.21iip + 00 6.450 £ ,07 0 OOOp.Ol 

2 ,2llp,00 0,45Oe*O7 0 ^OOO F .O1 

2«21lP + 00 6.450 E*o7 8,000E**01 


Rho gam t 

5,089E-°5 1,147 £+00 4,728E+°3 

5,O09 E -Os 1,147 e *00 4,72 be*03 

R, O09£-O5 1,147f*00 4,720E*O3 

b,°89E“^5 l,l47p,00 4,728£*03 

8,O09 E .O5 1 1 l47p + 00 4,72BE+°3 

t »»° 0 9E“®5 l,l47F+ f '^ 4,720£*O3 

S, 089E-°5 1.147F+00 4.728E+03 

b,0B9E-°5 1.147F+00 4,7p0E+O3 

S,069E-°5 1 , 1 47f + 0 0 4.720E + O3 

5,089E-°5 1.147F+00 4,720 E *O3 

5«0 69 £*05 1.147E+00 4«728E+o3 






I 

Z 



p 

Q 

PHE 


— l- 

0. 

-ul 

1 

200 r +00 

7.4?0g; + 02 

0. 

0 

-5". 00 °t« 

5 

000.. *02 

7 . 4?0 F «.0 3 

0 

• 

0 


o. 

.ul 

1 

?oo r *oo 

7 fc20p+°3 

O ' 

0 

2 

-4* 0 OOt M 

5 

0 0 Op* 02 

7,4?0 E *03 

o; 

0 

‘ 

0. 


1 

200^*00 

7 .420 *03 

0 

0 

3 

-3.000 


* 

0 0 01. * 0 2 

7* 4?0|*0 3 

0 * 

0 


0. 


3 

2 0 0 p * 0 U 

/,^20 £ *C3 

0* 

0 

4“ 

-2.000 

.01 

c; 

0 0 Op *02 

f 

0 

• 

0 


o. 


i 

?00p*00 

7.420 e *0 3 

0 

• 

0 

5 

-i.ooo t . 

.01 

5 

000p + 02 

7^420^+03 

0. 

0 

" 

0. 


1 

200 *00 

7*4?0 E+ 0 3 

0, 

0 

6 

0. 


B 

00(lp*02 

7.4?0g*0 3 

0. 

0 


0. 


1 

200 *00 

7.4 ? 0p+0 3 

0 

0 

7 

i;ooo t . 

.01 

B 

000p*02 

7#<f ?0 E *0 3 

o t 

0 


o. 


1 

200 *00 

7.420 E *0 3 

0, 

0 

0 

2,000 e . 

.01 

5 

0 0 Op* 02 

7.420 e +0 3 

0 * 

0 


0. 


1 

20 Op* 0 U 

7,420 E 4.0 3 

0 

• 

0 

9 

3,000 t . 

.01 

5 

000 F *02 

7.420 fc+ 0 3 

0. 

0 


0. 


3 

c o o p * o o 

7.420p*0 3 

0 

0 

. — ln 

4.000... 

.01 

c 

000p*02 

7.4?0 e *0 3 

0 * - 

0 

11 



l 

200 *00 

7.4-?0p + 0 3 

0 

1 

0 

5,O00 t „Ul 

5 

000p*02 

7,420 E+ 03 

0, 

0 
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M 

H 

phi 

2,21 lp*00 

6*45 0 £40 7 

e.ooOp-oi 

2,2llp + 00 

6.45°E*0 7 

8, OOOjp.Ol 

2,2ll F+ 00 

6,45°£ + 0 7 

8% 000 p .01 

2,21 lp+OO 

6,45°e*0 7 

0.000^01 

2,211 f *00 

6.450 e+ 0 7 

e# 000 E -01 

2,21 lp. + 00 

6,45 0 £+0 7 

8# 000p.01 

2,21 lp*00 

6.45°E40 7 

B.OOOp.Ol 

2,21 lp*O0 

6,450e+0 7 

8 ,ooo e -oi 

2,21 lp+00 

6.450£ + 0 7 

S,000 E -01 

2,21 lp*00 

6,45°e*0 7 

b # ooo e .oi 

2, 21 1 F + 00 

6,450 e+ 0 7 

s,ooo e -oi 
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rho 

GAM 

T 

5,0 8 9E-05 

1 ,l4 7 (T400 

4.728E*°3 

5,089E“°5 

1.147E* 00 

4.728E*°3 

5,°89e-°5 

i # U 7e;+ oo 

4 .72 8E *°3 

5, 0 89E-°5 

1,147e* 00 

4,7 2 8E+°3 

5,089e-°5 

1,147^00 

4 , 7 2 8 E ♦ ® 3 

5,°89e-°5 

1,147 E *00 

4,72 8e+ 03 

5,0 89 e-°5 

1 , l47 E *00 

4.72 8 E + f | 3 

5,089E-°5 

1, 147E + 00 

4.728E*03 

5.089E-°5 

1,147E* 00 

4,728E*03 

5,0 8 9E-°5 

1,147e* 00 

4»72 8 e+03 

5,089E-°5 

1 , 147e* 00 

4,7 2 8E* 03 



TR 166 - Vol. I 


Page 140 


COWL AND EPEE STREAM DA I A AT KOUNI - 70 


TABLE lib 


J = 1 

1 z P 0 PH6 St — w 


1 

-5.0006-01 

1 .639E + 02 

3.6806*03 

0. 

0. 

2.7716*00 

— 2 

-2.8166-01 

1.6376+02 

8.6816*03 

- 4.2416-03- 

-or — 

2. 7726+00 

3 

-5.9896-02 

1 .6306 + 02 

8.6846+03 

1 .029E-02 

0. 

2.7736*00 

^ 

” T. 63 66- Oh- 

1 . 6236*02” 

~“8.688E*03 

- 1.9536-02 

Oi 

2.775E+00 

5 

3.8866-01 

1.6166*02 

8.6936*03 

3.2646-02 

o. 

2.7786*00 

6 

6TTT3E -0 1' 

1*61 Op +“02" 

- 697F + 03 

4 * 9 7 7tr-Uc 

0 * 

"TTT B 0 tr^trO — 

7 

8.4846-01 

1 .605E + 02 

8.6996+03 

7.0576-02 

0. 

2.781E+00 

~ 8 

1 .0836*00' 

1.6006*0? 

* 70 3F +03 

9.44 1 E-02 

0 • 

Z « 783F + 00 

9 

1.3136+00 

1 .594E + 02 

8.7076*03 

1 .1986-01 

0. 

2.7646+00 

10 

1 .7786+00 

1 . 578n + 02 

3.71 76 + 03 

1 • 7456-01“ 

0 • 

2.7876*00 

11 

1 .7786 + 00 

1 .3656*02 

8.854E+03 

2.1596-01 

0. 

2.657E+00 

"1 T 

1 . 7786*00 

I • I S 3 F +13? 

9T009F+03 

2*o316 — 01 

0 . 

~2“r9'396“*~0 0 — 

13 

1 .778E+00 

9.3996*01 

9.1866+03 

3.1846-01 

0. 

3.0376*00 

T4 

1.7786*00 

9.3996 + 0 1 - -9. 1886 +03 

3.1846-01“ 

-~0. 

3.0376*00 

15 

1 .778E + 00 

9.3996+01 

9.0406+03 

3.184E-01 

0. 

4.1206*00 

1 6 

1 • 7786 +00 

9 . 399E +TJ 1 

9.0406*03 

3,1 84£-0r 

o . -- - 

4 . 120E + 00 

17 

1 .7786 + 00 

1 .0006*01 

9.9006+03 

0. 

0. 

6.7556*00 

18 

2.5006+0 1 

1 * 000E + 0 1 

"'"9.90 0E + 03 

0 . 

0 • 

6.7556* v 0 — 

19 

5.0006+01 

1 .000E + 01 

9.900E+03 

0. 

0. 

6.7556+00 


J-s— 2 ■ 

I Z P 0 

1 -5. 0 00E-01 1.6766+02 8.6576+03 frv 


PHE 


SI 

^ 8 - 7 - 0f - O 5- 

-3.3576-04 

■2.1426-04 

-3.0616-04 

■1 *9586-04 

-3.5066-04 

"p n C i\ /■ 

W C ♦ C 

-1.2146-04 
-3.6546-04 
•1 .1016-03 
►I .0646-03 
■ 1 .0656-03 
1.0456 - 0 3- 
-1.0456-03 

0. 

0. 

-<H- 


M 


-2t-7-60E^"00- 

2.7606*00 

2.764E+00 

2.7676+00 

2.771E+00 

2.776E+00 

2*76 0 6+03- 


2 -2 
— 3--S 
A 1 
5 — 3 
6 6 

—f 8- 

8 1 

-9 1 

10 1 

-H 1 

12 1 

-T-3 h 

14 1 

is — i- 
16 1 

17 - 1 

18 2 
-f9 — S 


.776E-01 
.41 2E-02 
.696E-01 
.9SSE-01- 
.2426-01 

• 063E+00 
• 315E*00 
. 778E *00 
.7786*00- 
• 7786 + 00 


1.6736*02 
1 .6646*02 
1 .6516*02 
1 .6376+02 
1 .6236*02 
1.6096 * 02 
1 .5966+02 
1 .5846*02 
1.5626*02 
1. 353E + 02 - 
1 • 143E+02 

9 . 3416+01 
' 9 . 3416*01 
9 . 3416*01 
1 . 0006 * 01 - 
1 . 0 o 0 1 ♦ 0 1 
1 . 0006*01 


8.6596*03 
8 ,6646*03 
8.6726*03 
-8.6806*03 
8. 6896*03 
8 .6 9 76 * 03 
8.7056*03 
8.7136*03 
8.726E+03 
8.8626*03 
9.0166*03 

9.1936*03 
9.0466*03 
9. 046E*03 
9.9006*03 
9.9006*03 
• 9.9006*0 5- 


6.5616-03 
hi 45 36 “02 
2.5036-02 
3.8736-02 
5.5846-02 
-7- . -59 0 6— 02 


.7786*00 
.7786*00 
.7786*00 
.7786*00 
.500E *0 1 
rft006*"0l~ 


9.861E-02 
-1.2326-01- 
1.7456-01 
2 .1566-01 
2.624E-01 
3. 1716 - 0 1- 
3.171E-01 
3.1 7 16— 0h 
3.171E-01 

0.- 

0. 

-9-i 


0. 
0 . 


2.784E+00 
2.7876*00 
2.7926*00 
2.862E*00 
2.9436*00 
-- 3 - yO 4O6*-00- 
3.0406*00 
4.1296*00 
4.1296*00 
6.7556*00 
6.7556*00 
-6 .7 55 6 *0-0- 
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5.435E*07 

5.4356*07 

5.432E*07 

5.4276*07 

5.423E*07 

b- . - 42 - i e ^7 - 

5.418E*07 

5.41 BE* 07 

5.4126*07 

5.404E + 07 ' 

5 . 283E+07 


8.000E-01 

8.000E-01 

8.000E-01 

8.000E-01 

8.000E-01 

8.00QE - 01 

8.000E-01 

8.000E-01 

8.000E-01 

8v000E-01 

8.0006-01 


S.OOOE-Ol 

8.000E-01 

0. 

0-. 

0. 


1.933E-05 
1.931E-05 
1.924E-05 
1.917E-05 
1 . 909E-05 

1 A C 

1.897E-05 
1.893E-05 
1.886E-05 
1 .870E-05 
1 .651E-05 

i i ^ a r 


1.197E-05 
1 . 197E-05 
2.565E-05 
2*565E-05 
6.456E-06 


1.1S9E*00 
1 • 159E*00 
1 • 1596*00 
1 • 159E+00 
1 • 1596*00 
-ii459E*-O0 
1 • 159E*00 
1 • 1596+00 
1.1596*00 
1 .1596*00 
1 • 161E+00 

1 1 r . A A 


1 . 166E+00 
1 • 166E+00 
1.3146*00 
1 .314E+00 
1.387E+00 


9.9506*06 


6.456E-06 1,3876*00 


4.218E+03 
4.217E+03 
4 .2 16E + 03 
4.214E+03 
4.212E+03 
-4.21 1E + 03- 
4.209E+03 
4.208E+03 
4.206E+03 
4 . 202E+03 
4.145E+03 
-0796+03- 
4.0046*03 
4. 004E + 03 
2. 1356*03 
2.1 35E+03 
9, 024E+02 
-9t-024E-*02- 
9 . 024E+02 


5.453E+Q7 

5.4496+07 

5.442E+07 

S.435E+07 

5.427E+07 

5 .4 1 3E* 07 
5 • 4 0 7 E ♦ 07- 
S.395E+07 
5.275E+07 
5* 1 J8E + 07 
-4- .977E + 07 
4.9776*07 
1 .8046+07- 
1 .8046+07 
9 .950E+06 
9.950E+06 


8.000E-01 
8*0006-01 
8.000E-01 
8 • OOOE-0 1- 
8.000E-01 

-8-W00- O 6— 01- 

8.000E-01 
8*0006—01- 
8*000E-0 1 
8. 00 06-01 
8 • OOOE-0 1 

a r \ rv n c- — n) 

“ IUVUC. v * 

3.000E-01 
0 . 

0* 

0* 

0. 


1.967E-05 
1 . 959E-0S 
1.945E-05 
1.931E-05 
1.917E-05 
4^9046-05- 
1.890E-05 
4.877E-05 
1.854E-05 
- 1 *6386-05 
1 *4176-05 


1.1916-05 

2,5586-05 

2.558E-05 

6.456E-06 

6.456E-06 


1 • 159E+0 0 
1 • 159E+00 
1 • 159E+00 
1.159E+00 
1 • 159E+00 
4-.459E* 00 
1 * 159E+00 
1 • 159E+00 
1.159E+00 
1.1616*00 
1.1636*00 
4-^1866+00 
1 *166E*00 
1,3146*00 
1.3146*00 
1.3876*00 
1.3876*00 


-4.-2 276 + 03- 
4.2266*03 
4.224E+03 
4.221E*03 
4.217E*03 
4.214E+03 
4.218E+03- 
4.207E+03 
4.2Q4E+03 
4.198E*03 
4.1416*03 
4.076E+03 
-4.002E+03- 
4.0026*03 
2. 128E+03 
2. 128E*03 
9.024E+02 
9.024E*02 
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J = 3 

1 Z P 0 PHE S-I — M 

1 -5.000E-01 1 . 667E+02 8.663E+03 0. -1.329E-04 2.762E*00 

2 ~~ -2. 827E-0 1 1 ,664E*02 8 7 6646 ♦ 03 4v42fr t r -~ 9 3-- » . 6 44E - 0 5 2 -g- 7 6 4£- »- 00 

3 -6 . 349E-02 1.65SE+02 8.670F+03 1.059E-02 -2.453E-05 2.766E+00 

4 1T615E-01 1.643E>02 — ftvfr77E+03 — 1 .-985E-02 — S.974E-05 2.770E+00 

5 3.811E-01 1.630E+02 8.684E+03 3.245E-02 1.260E-04 2.774E+00 

6 6~. 1 176-01 W61 8E ♦02- 8.6916^03 — 4.94 3E -02 9.871 E-OS 2 .777E+00 

7 8.403E-01 1.607E+02 8.698E*03 6.960E-02 6.624E-05 2.761E*00 

- — 8 1 ♦ 0 7 1 E ♦ 0 0 i .599E * 02 8.703E * 03 9- . - 278E - 02 5 . - 1 29 E - 0 5 2 .7 H 2F + 0 9 

9 1.306E+00 1.S91E+02 8.707E*03 1.186E-01 1.639E-0S 2.785E*00 

— 10 “1 .778E + 00 1. 576£ ♦ 02 8Vtt 76+03 H-74SE-01 — 4 . 14 IE -OS 2 . 788E + 0 0 

11 1.778E*00 1.364E+02 8.854E^03 2.15BE-01 -4.077E-05 2.8S8E*00 

12 K77BE+G0 — K152E^02— 9.009E + 03 -2.6306-01- -4^.00 7E -05 2.940E+00 

13 1.778E+00 9.393E+01 9.187E+03 3.182E-01 -3.929E-05 3.037E*00 


H 


-PH6 


-RHO- 


-GAM 


S.4S0E+07 8.000E-01 1.961E-05 1.159E+00 


S.444E+07 
5,438E*07 
S.431E+07 
5.425E+07 
5.420E+07 
5. 4 15E + 07 
5.41 OE+07 
-5-.403E+07 
S.283E+07 
S,144E*07 
4.982E+07 


8.000E-01 
8.0006-01 
8.000E-01 
8.000E-01- 
8.000E-01 
-frrOO OE - Ol 
8.000E-01 
8.000E-01 
8.000E-01 
8.000E-01- 
8. Q00E-01 


1.949E-05 
1 .9376—05 
1.924E-05 
-1.912E-05 
1.901E-0S 
1 _ 89 1 F- Q 5 
1.884E-05 
1.870E-05 
1.650E-05 
1.426E-05 
1.197E-05 


1 . 159E+00 
1 • 159E+00 
1 • 159E+00 
-1.1 59E + 00 
1 • 159E*00 
E-+00- 
1 • 159E+00 
-1.1S9E*00 
1 • 161E+00 
-1 • 163E*00 
1 • 166E+00 


4 . 225E>03 
>03- 
4 .222E + 03 
4.219E+03 
4 . 2 1 6E*03 
4 .2 1 3E+03 
4.21 OE+03 
-9h.208E-»O3- 
4.206E+03 
4 .202E* 03 
4. 144E*03 
4.079E+03 
4. 004E*03 
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14 

1. f fttfc*U0 

9.3S>Jfc>Ui 

9. Ib/fc+UJ 

J. 


J.UJft+WW 

15 

1.7786*00 

9.3936*01 

9.0416*03 

3.1826-01 

0. 

4.1216*00 

“ 16 

“T. 7786*00 

9.3936*01 

~Q.04TE*03“~ 

3 71826-01 

0 . 

4.1216*00 

17 

1. 7786+00 

1.0006*01 

9.9006*03 

0. 

0. 

6.755E+00 

IP 

2.5006*01 

' 1.0006*01 

9.9006*03 

0. 

0 » 

b • 755E ♦ 00 

19 

5.0006*01 

1.0006*01 

9.9006*03 

0. 

0. 

6.7556*00 


i 

J = 4 

2 

P 

0 

PHE 

SI 

M 

i 

-5.0006-31 

“1.6766+02 

8 » 657E *03 

0. 

t> « 7 bBt -03 

2 • 760E+00 

2 

-2.7826-01 

1.6736*02 

8.6596*03 

6.4436-03 

4. 5026-03 

2.761E+00 

3 

-5T.48TE^0T 

“TT^646TiJ2- 

8 • 664E ♦ 0 J 

1 .4 316-02' 

4.0^96-03 

Z » 7b4t *00 

4 

1.6966-01 

1.6516*02 

8.6726*03 

2.4846-02 

2.7456-03 

2.7686*00 

5 

“3.9506-31 

“T. 6376* 02 

“8". 68 06* 03 

“3.853F-0? 

"3.7986-03 

2.7 72E*00 

6 

6.2436-01 

1.6236+02 

3.6906*03 

5.5786-02 

9.1036-04 

2.7756*00 

7 

8753 06 -Or 

“1.6 096*0 2 

8 « 697E + 03“~ 

“7T5976-02 

-1-.-99SE-04 

2. 779E+00 

8 

1 .083E + 00 

1 .5966*02 

8.7056*03 

9 .8636-02 

-7.6816-04 

2.7846*00 

9“ 

1T3I5F*W~ 

1 .5846 + 02 

9 ".71 36*03 

1 . 232t-0 1 

-2.479E-03 

z • 

10 

1.7786*00 

1 .562E + 02 

8.7256*03 

1 . 7456-01 

-8.864E-03 

2. 7926*00 

11 

1.7786*00 

173536*02 

“81861 E+ 03“ 

273566-0 I" 

•-8 77286-03 

2.862E + 00 " 

12 

1.7786*00 

1.1436+02 

9.0156*03 

2.6246-01 

-8.579E-03 

2.9426*00 

o 

1.7786*00 

9.342 6 ♦ 01 ' 

“971926*03 

3.1716-01 

-8.4146-03 

3.0396*00 

14 

1.7786*00 

9.3426*01 

9.1926+03 

3.171E-01 

-8.4146-03 

3.0396*00 

nr 


■L-mnaain 

■■•HiLi.mOM 

3.1716-01 

0. 

4 ;■ 1 

16 

1.7786*00 

9.3426*01 

9.0466+03 

3.1716-01 

0. 

4.1296*00 

17 

“1T778E + 00 

1 .OOOE + 01 

‘ 9.-9006*03“ 

“07 

0 • 

b « rbbr ♦ UU 

18 

2.500E*01 

1 .0006 + 01 

9.900E+03 

0. 

0. 

A- 

6.7556*00 

" 19 

“5.0006*01 

1 .0006*01 

- 9.9006+03 

0 * 

0* 

o * not ♦ uu 


J = 5 

— 1 2 P ' 0 PH£ SI M 

1 o. 1.5596+02 8.6496*03 -3.6546-03 -1. 4126-01 2.7956*00 

2 — 3.0006-01 — 1*; 5 74 6*02 — R. 6386*03 — 1V516E-14 -1.5256-01 2.7966*00 

3 3 . OOOE-O 1 1.2626*02 R. 8476*03 6.3116-02 -1.4906-01 2.902E+00 

— 4 3V0 006^03 9.50 76 *01 9 .10 0 6* 03 1.4116-01 -1.4496-01 3tO386+-&0 

5 3. 0006-01 6.3936*01 9. 4286*03 2. 4496-01 -1.3996-01 3.2276+00 

— 6 3.0006-OT 6 i 3936*01- — 974286 *03 2 . 4 4 96 -0 1 — 3 v 3 996 -0 1 3.2276*00 


7 3 . OOOE-O 1 6.3936*01 9.335E+03 2.448E-01 0. 4.659E+00 

“8 — 370006-01 "6.39 36* 01 9 . 3356*03 2.4486-01 0; 4.6596*00 

9 3. 000E-01 1.000E+Q1 9.9006*03 0. 0. 6.755E+00 

rtr — 2.50 06*01 1 .0006 + 01 9.9006 + 03 07 07 fcrv?55E*00 

11 5.000E+01 1.0006*01 9.9006*03 0. 0. 6. 7556+00 


1 

2 
3” 

4 

5 

6 

~n~ 

8 

9 

10 

rr 


J - — 6~ 

1 

751 

3.0006-01 
“310 00 6-0 T" 
3.000E-01 
37 0006-0 1" 
3.0006-01 
“370-006=^0- 
3.0006-01 
"T. 0006-01 
2.5006+01 
5.0006+01 


1.5936+"02“ 
1.5926*02 
07 2766+02 
9.5916+01 
-674266 + 0 1 
6.4266+01 
b. 4266* 01 
6.4266*01 
1.0006*01 
1.0006*01 
1.0006+01 


0 PHE 

R . 69 1 6*03 3.049 6 

4.6836*03 -1.0066 
""87 8936*03 — 673216 
9.1476*03 1.414F 

9.4766*03 ' '27458E 
9.4766+03 2.4586 

-97331T*03 2.4576 


M 


SI 

• 0 4 -7.3006-02 2-7735 6 * 00 - 

•14 -8.9746-02 2 . 786E+00 

-02 - 8.7 636- 02 — 2.8936*00 
-01 -8.5216-02 3.0306*00 

-01 -8.2266-02 - 3.2216*00 
-01 -8.226E-02 3.2216+00 

4.6516 *00-- 


9.331E+03 2.4576 

979006*33 0 . 

9.9006*03 0. 

9.9006*03 Ot^ 


=03- 

-01 


0. 

-trr 

0. 

o-i- 


4.6516*00 
6.7556+00 
6.7556*00 
4. 7SRf*00 
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l.B09E*07 0, 2.564E-05 1.314E+00 2.134E*03 

J.B09E+07 — Oi “2“i564E~05 1.314E*00 2.134E+03 

9.950E*06 0. 6.456E-06 1.387E+00 9.024E+02 

- 9.950E + 06 '0¥ 6-.456E-06 - 1.387E + 00 9.024E*02 

9.950E*06 0. 6.456E-06 1.387E+00 9.024E*02 


H 

S.4S4E+07 - 

S.453E*07 


PHI 

8.UOOE-Or 

8.000E-01 


RHO 

1.970E-05 

1.967E-05 


GAM 

1 . 159E + 00 
1 . 159E + 00 


4.227E*03 

4.226E*03 

-4t224E>03- 


5.442E'M)7 B.OOOE-Ol 1.945E-05 1.159E*00 4.221E*03 

-_ 5'.435E* 07 — »;OOOE-Ol 1.931E-05- 1.159E+00 4.217E*03 

5.427E+07 8.000E-01 1.916E-05 1.159E+00 4.214E*03 

"Sv4 19E ♦ 07 — 8T000E-0 1 — 1 • 902E-05 — 1 . ! 59E ♦ 0 0 4 . 2 1 0 E *03 

5.4l3E*07 8.000E-01 1.890E-05 1.159E+00 4.207E*03 

5.407 E +07 8»000£*01 1.877E-0S lrlS9e*90 4i204E*03— 

5. 395E+07 8.000E-01 1.854E-05 1.159E+00 4.198E*03 


“ 5t 275E*07 — 8.000E-0 1 lv638E**05 " 1 • 161E + 00 4.141E+03 

5.1 38E*07 8.000E-01 1.417E-0S 1.163E + 00 4.076E*03 


-4^977E*97 — 8.000E-01 I . 191E-05 1 . 166E + 00 4.002E+03 

4.977E*07 8.000E-01 1.191E-05 1.166E+00 4.0l)2E*03 

1.804E«07 - - - Ow 2.55 8E~95 1 4E+00 — 2 .12 8 E+ 0 3 

1 . 804E+07 0* 2.558E-05 1.314E+00 2.128E+03 

-9i950E*06 ~ 0. 6.456E-06 1.3B7E+00 9.024E+02 

9.950E+06 0. 6.456E-06 1.387E+00 9.024E+02 

9.950E + 06 — O* G.456E-06 1.387E+00 9.024E*02 


-H- 


PHI- 


5.388E+07 8.000E-01 

-5r395E ♦ 0 7 — 8v000E-0 1- 
5.2 12E + 07 8.000E-01 

4;-9Q5E * 07 6. Q 00 E-01 

4.682E+07 8.000E-01 

‘♦.662E*07 — 8.000E-01 
1 ,539E*07 0. 

1 .539E+07 — 0. 

9. 950E+06 0. 

9.950E + 0 6 9-i 

9.950E+06 0. 


RHO 

1.851E-05 

1.867E-05 

1.543E-05 

1. 2 09E-05- 

8.606E-06 

8.606E— 06 

2. 138E-05 

2.438E-05 

6.456E-06 

6 * 45 6E - ( K>- 


6.456E-06 


GAM T 

1.159E+00 4.196E+03 

1 . 159E+00 4. 199E+03 

1.1 62E+00 4 , 1 1 2E*03 

;-H66E+00 — 4-,006E^-08- 
1 . 1 72E+00 3.871E+03 

1 • 1 72E+00 3.871E*03 

1 « 343E + 00 1 • 742E* 03 

1.343E+00 1 . 742E+03 

1 .387E+00 9. Q24E+02 

1.387 E+ 0 0 9-.-024E*02- 

1.387E+00 9. 024E+02 


H PHI RHO GAM T 


5.406E*07 

5t222E*07- 

4.993E*07 

4.686E+07 

4.686E+07 


8.000E-01 

8»000E“0l 

8.000E-01 

-8.000E-01 

8.000E-01 


1 • 885E-05 
1.557E-05 
1.218E-05 
8.638E— 06 
8.638E-06 


1.159E*00 
1 . 162E + 00 
1 . 166E + 00 
1 * 1 7 1E+00 
1.1 71E+00 


4.204E*03 
4. 1 16E*03 
4.010E+03 
3.873E*03 
3.873E*03 
-l-.*746E-*-03 


1 ,542E*07 0. 2.144E-05 1.343E«-00 I.746E*03 

~9.950E*06 — 0. 6.456E— 06 • 1.387E+00 9.024E + 02 

9.950E+06 0. 6.456E-06 1.387E*00 9.024E*02 

9.950E*06 0. 6.456E-06 -1.387E+00 9.024E*02 
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J = 7 


I — 

r 

P 

- ■— tJ- 

— PHE —Si 

M 

1 

0 . 1 

« 633E ♦ 0? 8 

.6876*03 2 

.305E-03 -2.4386-02 2 

.7736+00 

~z~ 

3 . 000 F -01 - 1 

V63?E^02 — 3 

.692E*th3-~2 

•i 77E-14— -3.-47 16-02' 2 

.7766+00 

3 

3.000E-01 1 

.305E>02 8 

.9036*03 6 

.371E-02 -3.388F-02 2 

. 8836*00 


H 1 

PHI 

RHO 

GAM T 



5.4276*07 

8 . 000 E -01 

1.927E-05 

1 . 1596*00 4.2146*03 



5.4276*07 

8 . 000 E -01 

1.926E-05 

1 . 159E*00 4,2 146. + 03 



5.24 IE* 07 

8 . 000 E -01 

1.589E-05 

1 . 162E*00 4.1 256+03 


u 

3.000E-01 

9.788E+01 

9.1606*03 


3.0216*00 

5 

3.000E-01 

6.522E*01 

9.494E+03 

2.4R4F-01 -3.178E-02 

3.215E*00 

6 

3.000E-01 

6.522E*01 

9.4946*03 

2.484E-01 -3.178F-02 

3*215E*00 

7 

3.000E-01 

6.522E*0 1 

9 . 322F *03 

2.483E-01 0. 

4.631E*00 

~8 

3.000E-01 

6.522E*01 

9.3226+03 

2 .483E-0 1 - Q-- 

4.631E*00 

9 

3.000E-01 

1 ,000E*01 

9.900E+03 

0. 0. 

6. 7556*00 

nr 

2 .5OOE + 0 1 

“TTO'CTOE+0 1 

~9T90FOE + O3 

0 * 0 • 

5.TS5E*0(T 

n 

5 » OOOE + 0 1 

1 ,000E*01 

9.900E+03 

0. 0. 

6.755F+00 




- 




J = 8 





i 

z 

P 

Q 

PHE SI 

M 

i 

0 . 

■ ITOMMIM 

8.6766*03 

5 ■ 768E-03 0"J 


2 

3.000E-01 

1.663E*02 

8.6726 *03 

8.390E-1S 0. 

2.7646+00 

3 

* “3.000E-01 

1. 3296*0? 

"8.88 6F>"03" 

6V420?-02 O - . 

2.872E+00 

4 

3.000E-01 

9.948E+0 1 

9. 1456*03 

1.439F-01 0. 

3.01 IE* 00 

■ - 5 - 

-tjtoooft-ot" 

6.6076*01 

9.4836*03 ' 

2.507E-01- 0. 

3.2066*00 

6 

3.000E-01 

6.607E*0l 

9.4836*03 

2.507F-01 0. 

3.206E*00 

/ 


tJ*f>07t£ + 01 

9.314E+03 

2.506E-01 (T; 

'4 . ni3E ^01 r 

8 

3.000E-01 , 

6 . 607E + 0 1 

9.3146*03 

2.506E-01 0. 

4.613E+00 

9" 

3.000E-01 

1 • 000E> 0 l 

9.9006+03— 

0 . 0 

b • 755E + 00 

10 

2. 500E *0 1 

1 .0006*01 

9.9006*03 

0. 0. 

6.755E*00 

— n 

5. OOOE *0 1 

1 .0 0 061 + 0 1 

9.9006*03- 

Ov 0. 

6.755E*00 


5.009t*0 t 

a^uuut-ui 

1 .24.ue.-us 

l . lose. ♦ uu " 



4.697E+07 

b.oooe-oi 

8.753E-06 

1 . 171E+00 3.878E*03 



4.6976*07 " 

8.000E-01 

8.753E-06 - 

1 • 1 7 1E+00 3.878E+03 



1 .551E*07 

0. 

2* 160E-05 

1 .342E *00 1 .759E+03 



1 .581E*Uf 

0* 

'2*1606-05 

1.342E+00 1.759E+03 



9.950E*06 

0. 

6.456E-06 

1 .3876*00 9, 0246*02 



9 *9d0E*0o 

0 « 

6 . 456E-06“ 

1t3876 + 00— 9. 0246*02 — 



9.950£*06 

0. 

6.456E-06 

1 .3876+00 9. 024E+02 



H PHI RHO GAM T 


r.-159£ *00 — 4 i 220E+03-- 


S.443E+07 0 ■ OOOfc'-0 1 1.9S7E-05 1.159E+00 4.222E+03 

5v255E*07 “8.000E-01 1.612E-05 1.162E*00 4.132E*03 

S.021E*07 8.000E-01 1.257E-05 1.165E*00 4.023E+03 

4.706E+07 -8.000E-01 8.845E-06 1.1716*00 3.882E*03 

4. 706E+07 8.0006-01 8.845E-06 1.1716*00 3.882E*03 

1 .S SSfrrftt — (hr- 2-r4-74t-05 Fr34i€*<Kl H'77 1E*03 

1 .558E+07 0. 2 • 1 74E-05 1.341E*00 1.771E*03 

9.9506*06 0* 6.456E-06 1.3876*00 9.024E*02 

9.9506*06 0. 6.456E-06 1.3876*00 9.024E*02 

9.9506*06 — 0* 6.456E-06 1.3876*00 9.0246*02 
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ALPHA 



1 4.347E-01 0. 3.184E-01 0. 0, 0. 0. 


' —2 4v334E--0i Ov- 3^ 171 E -01-0^ HK 0. 0. 

3 4.346E-01 0. 3.182E-01 0. 0. 0. 0, 


2.449E-01 0. 

^.458E**01 (hr 

2.484E-01 0. 

“2t507E“01 0~. 
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KOUNI - 

120 







TARI F Tic. 




Z MOMENT 

AXIS = -0. 

X MOMENT AXIS = 

-0. 

PITCHING 

PITCHING 

MOMENT *" 
MOMENT = 

“ l.I24E*04 
1 • I 05E*04 

THRUST" '= 
THRUST = 

37 810E*02 

3.8106*02 

“LIFT"- 1.3466*03 

LIFT - 1.3336*03 

THRUST - 

3.81 06*82 

LIFT = ’ 1 • 3466*03 

PI ICHlNO 

MOMENT — 

m24E-*04 



J = 1 



Y = 

0 . 


I 

z 

p 

Q 

PHE 


M 

1 

- 5 . OOOE-Ol 

2 . 240 E *01 

1 . 0206*04 

0 . 

0 . 

3 . 7246*00 

2 

- 1 . 3816-01 

2 . 2326*01 

1 . 0216*04 

' 279 ' 18 E-T) 2 - 

0 . 

3 . 7256*00 

3 

2 . 3596-01 

2 . 2176*01 

1 , 0216+04 

5 . 846 E -02 

0 . 

3 . 728 E *00 


6 . 2476-01 

2 . 1956*01 

I . 0226+04 

7 ) . 7036-02 

0 . 

3 • 733 E +0 0 

5 

1 . 0296*00 

2 . 179 E *0 1 

1. 0226 + 04 

1 . 1306-01 

0 . 

3 . 737 E *00 

6 

1 . 4 - 5-3 6 ^ 00 “ 

2 . 19 OE + 01 

r: 75 ‘ 22 mpr 

1 , 345 E- 0 i 

U • 

“ 377356 * 0 0 “~ 

7 

1 . 8946*00 

2 . 2536*01 

1 . 020 E + 04 

1 . 495 E -01 

0 . 

3 . 7216*00 

a 

2 . 3456*00 

2 . 3976 * 01 ' 


0 . 

3 *o 91 E+U 0 

9 

2 . 7836*00 

2 . 6306*01 

1 . 0106*04 

1 . 558 E -01 

0 . 

3 . 6446*00 

— ro- 

“ 3 . 1816*00 

27945 E* 0 T 

1 . 0036*04 

175746 - 0 - 1 “ 

0 • 


il 

3 . 597 E *00 

3 . 336 E +0 1 

9 . 9396*03 

1 , 6756-01 

0 . 

3 . 526 E +00 • 

rpr 

3 . 8146 * 00 " 

“ 3 . 5506*01 

9 . 8976*03 


0 . 

3 . 495 t*UU 

13 

3 . 814 E *00 

3 . 5506+01 

9 . 3856*03 

1 . 762 E -01 

0 . 

4 . 7696*00 

14 ~ 

4 , 2516*00 

4 . 4736 * 0 1 " 

' 9 . 3476*03 

2 • 0566 - 01 “ 

0 • 

4 . 6806*00 

15 

5 . 0546*00 

6 . 421 E *0 1 

9 . 3326*03 

2 . 4566-01 

0 . 

4 . 6536*00 

16 

5 . 0546*00 

1 . 0006*01 

— 979006+03 

- 477676 - 14 - 

07 

677556*00 1 

17 

2 . 5006*01 

1 . 0006*01 

9 . 9006*03 

0 . 

0 . 

6 . 755 E +00 ' 

nr 

“ 57 cnT 0 E'*Trr 

1 . 0006 * 01 “ 

9.9006 + 03 " 

Oa 

0 * 

6 • r boE + tru 


t~ 

1 

sr 

3 

“4“ 

5 

~6 


-5. 

^r. 

2. 

6 . 

1. 

r. 

1. 


J = 2 

— Z 

000E-01 

rsGE^trr 

7376-01 

6906-01" 

0776*00 

4996+00 

9246*00 

3526*00 

7796*00 

1656*00 

5696*00 

7876+00 

787E*00 


2. 31 1 6 + 0 1 
2.3046+0T 
2.277E+01 
ZW 246ET+01 
2.226E+0 1 
2.2356*01“ 
2.291E+01 
2.A22ET^r 
2.6466+01 
2. 9456*01 
3.3156*01 
3.5266*01 
3.5266*01 
4.4266+01 
6.3476*01 
■IvOtJOE^Ot- 
1.Q00E+01 
1.000E*t)l- 


XJ 

1 • 01 7E+04 
1 .0176 + 04 
1.0186*04 
IV019E+04— 
1 .020E+04 
"17 0206*04 

1 .018E + 04 
i;e rr5E^TTA— 
1.0096*04 
1.0026+04“ 
9.9406+03 
9V899E+03 _ 
9.3896 *03 
r*93- 
9.3386*03 
“979006*03-- 
9.9006*03 
9.9006+03“- 


PH 6“ 

0. 

3. -2 T6 e 


3.000006-01 

SI 

3.7606-02 
— 3T-7t^f^2— 


6.339E- 
972156' 
1.1786' 
1V385E- 
1.5146' 
1 »557 E~ ' 
1.5436 
l t554F- 
1.6476 
1.7306 
1 . 7356 
2.027E 
2.4316 
-4t6026~ 
0. 

0 . 


= 02 " 

-02 3.6736-02 

-02- 3.5876-02 
-01 3.S19E-02 

•Ot — 3 ,4666-02 
-01 3.492E-02 

^tn — 3-.- s -7t)r -0 2- 


M 

3.7 13E+00 
“ 3.7146*0 0 — 
3 .7216*00 
3.7286*00 ~ 
3.7326*00 
3.7326*00 
3.71 9E* 00 
3.688 6^00 — 
3.640E+00 
3.5876*00 - 
3.5296*00 
3.4986*00 
4.7806+00 
4.098E'H)0 


“8 27 

9 2. 

10 — 3. 

11 3. 

12 3. 

13 3. 
14 — 47 


15 

"re- 

17 

18 


5. 

“57 

2 . 

5. 


032E *00 
0326*00" 
500E+01 
OO 06*01 


■01 3. 5866-02 
■01 - 3.0866-02 
■01 2.5746-02 

■01 -2.0716-02 
■01 1.0456-02 

-04 1.3846 - 82 - 

■ 01 1.4966-02 
- H*—* 1 .4396- lb 
0. 


4.6696 4-00 
6. 7556*00 
6.7556*00 
6.7556*00 


O 
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R 

3.724E+00 
^J7725E*00 
3.728E+00 
3T733E«-O0 
3.737E*00 
3.735E* 00“ 
3.721E*00 

—3^691 E>0tr 

3.644E+00 

3i"587E*00 

3.526E*00 

3.49SE+00 

4.769E+00 

-4t68OE*O0“ 

4.653E+00 

~6t755E>00 

6.755E+00 

6.755E+00 


H 

3.996E*07 

3.994E*07 


pm 

8.000E-01 

-8T000E-01' 


3.990E+07 8.000E-01 
3.983E+07 — 8vOOOE=Or 
3.978E+07 8.000E-01 


3.983 F » TT7 8. 0 0 0E - 0 1 


4.001E+07 

4.041E*07 


B.OOOE-Ol 

8.000E-01- 

8.000E-01 


4, 101E+07 
4vl75E«-07 — 8iOOOE-Ol 
4.263E+07 8.000E-01 


1 *49 IE* 07 0. 

t^527E*07 — 0-.- 
1.S41E+07 0. 

9 ,950E*06 0. 

9«950E*06 0, 

- 9;9 5 fl£ - yft 6 — Ott 


RHO 

3.S29E-06 
3.519E-06 
3.500E-06 
-3.470E-06 
3.448E-06 
— 3T4 *>4F=06- 
3.546E-06 
3* 739E-06 
4.043E-06 
-4.449E-06 
4.948E-06 

j#tI tf tr~u o 

1.234E-05 
1.506E-0S 
2* 143E-05 
r.-456E-06 
6.4S6E-06 
6.4 5 6E - 0 G- 


GAW 

1.184E*00 
U184E>00 
1.184E+00 
1.184E+00 
1.184E*00 
"It 784 £“*^0“ 


1 • 183E+00 
1 * 182E + 00 
1.161E+00 
1 • 1 80E + 00 
1.179E*00 
-7v778E^OO- 
1 • 346E+00 
1.344E+00 
1 . 343E+00 
1 »387E*00 
1 .387E+00 


,599E*03 

♦598E*03 

,594E*03 

,587E»03 

,582E^03 

•'587E»03 _ 

*603E*03 

•617E*03 

,639E*03 

,665E*03 

.698E+03 

.714E*03 

,672E*03 

.724E*03 

• 746E + 03 
.024E+02 
. 024E*02 

• 024E'*02- 


M- 


-PRf- 


-PHO- 


3.713E+00 4 • 020E+07 8.000E-01 3.638E-06 


3.721E+00 
-3t 728E*00“ 
3 * 732E + 00 
-3v732E*00 
3.719E+00 
-3v 688E*00 
3.640£*00 
-3v587E*00 
3.529E+00 
-3y498E*00 
4, 780E+00 
4.6 9 8£ * (K h- 
4.669E+00 
6* 755E + 00 
6.755E*00 
-6*755E»00 


4.010E+07 8.000E-01 

-4 . 0 0 0 E ♦ 07 — 8 . 0 0 0 E- 0 1 
3.994E*07 8.000E-01 

3.997E*07 8.000E-01 

4.012E+07 8 • OOOE-O 1 

4. 104E*07 8.000E-01 

4vl 74E+07 — “B * 000E“0l _ 
4.2S9E+07 8.000E*01 

4.300E+07- 8.000E-01- 
1 .487E + 07 0* 


1 . 5 1 9 E + 07 9-»- 

1.535E*07 0, 

-9.950E+06— (hr 


3.592E-06 

3.550E-06 

3.524E-06 

3.537E-06 

3.609E-06 

3.777E- OG- 

4 • 064E-06 

4.448E— 06 

4.921E-06 

5.184E-06 

1 • 230E-05 

i C r>nr— n c 
1 4 5v Ut 


“GAM 

1 • 183E+00 
1. 163E*00 
1 • 183E+0Q 
1 .183E*00 
1 ,184E*00 
1. 183E+00 
1 ,183E»00 
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— 9.391 E*07 — 8r(XK>E-0 1 — 5 .806E-0 & trl7 6 E ♦ 0 0 ~ 3 • 74 8E ♦ 0 3 

— 1.417E ♦ 07 — 0* 1t506E— 05” 1 *352E*00 1.561E+03 

i>421£»07 — (hr 1.690E-Q5 — 1 ,35 IE* 00— H 5670 *03 

— 1t 549E ♦ 07 — 9* 2vi58E-05 1*3420*00 l-.-758E*03 - 

— 9.950 E*tf 6 — 0* ~~ 6.456E-06 1*387 E* 0 0 9 . 0 24E *02 

9 > 9 50E * 6 6 — fr* 6 . 4 56E - 0 6 1 . 3 8 70 *00 — 9*02 4 E * 0 2 

6t456E— 06 — 1*387 E ♦ 0 0 9. 024E*02 


— 9v950E'*06 — 0* 




4.108E-02 


6.243E-02 0. 

o. rtrr 


BETA 


1 3.577E-01 0. 

2 3.563E-01 07 

3 3.562E-01 0. 

4~“'3.5?BE-0r 0. 
5 3.180E-01 0. 

*6 372 ft 9E‘^n 0 . 

7 3.328E-01 0. 

8 ' 3.603E-01 07 


1 16 
2 — 16 
3 16 
A 16 
5 9 


7 9 

8 9 

3»483E*00 
6 • 397E-0 1 
3.782E-02 


0 

“ 0 — 

IIIII 

urn 

mu 


1.762E-01 

r.740E-01 

1.780E-01 

1.762E-01 

l.<*21E-01 



2.22SE+00 
■57528E-0r 
•9.898E-02 


IIIII 

IIIII 

IIIII 


1 1 1 1 1 
I I I II- 
IIIII 


0 

— a 

IIIII 

IIIII 

mu 


o 

o- 
0 

- 0 
0 

— t- 

o 
o 

mi 

nn 

HIM 
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EXTERNAL WRAP AROUND REGION 


I = 7 

R TH ~~ — P 0 PH E SI 

3.9S6F-01 0* 2.408F+01 1.009E+04 1.553E-01 1.017E-01 

3.8 7 2E-0 I 1 7 7455-01 ZTttVTTVl 1.008E+04 1 .51VE-01 771977-02 

3.787E-01 3.491E-01 2.473F*01 1.006E*0<* 1.485E-0I 4.293E-02 

3 » 7 0 3F - 0 1 5 • 2 36E - 01 2 . 5 0 SE ♦01 T • 00 5£ ♦ 0 4 '”"1 • 4 52E- 0 1 1.452E-02 

3.619E-01 6.981E-01 2.536F+01 1*003E*04 1.420E-01 -1.325E-02 

3.535E-01 8 . 7 2 7 E - 0 1 2 .5 6 6 F * 0 1 1.00 2 E ♦ 9 4 1.388F-01 -4.039E-02 
3.451E-01 1 • 047F + 00 2.596^*01 1.001E+04 1.357F-01 -6.688F-02 

3736 7E-01 rT727ET7TT5 — 27675 F'+TH 9T.993F + 03 1 . 32 7E-0r^T274t-l)2 

3.283E-0 1 1.396E+00 2.653F*01 9.980E*03 1.298E-01 -1.180E-01 

3.199E-01 ' 1 .571E+00 ~ 2.680E+01 9V968E4U3 ~I.270E-7jr -I7425E-01 


I 

= 8 





R 

“Th 

P 

u 

PHE 

51 

7.909E-01 

0. 

2.482F+01 

1 .007E+04 

1 . 569E-0 1 

1.082E-01 

7. 74 IE -01 

1.7456-01 

2.521 5 + 0 1 

1~. 0056* 04 

V.538E-0 1 

7.762E-02 

7.573E-01 

3.491E-01 

2.559E+01 

1.004E+04 

1 » 50PE-0 1 

4 . 7 70E-02 

7.40AE-01 

5.2366-01 

2.596F+0T 

1 .0026+04 

1.479E-01 

— 1 .844E-02 

7.2365-01 

6.981E-01 

2.6325*01 

1.001E+04 

1 .451E-01 

-1.016E-02 

7706ffE-7TT 

“5772 76=771” 

2766 /F + 0T”~ 

y .Rv4t + 03 

1 .423E-0T 

■"-3L8 r or- 0-2“ 

6.900E-01 

1 .0475+00 

2.7026*01 

9.979E*03 

1 .396E-01 

-6.539F-02 

6.7325-C 1 

1 .2226+00 

2.735E*01" 

9^966E >03 

1.369F-01 

-9.202E-02 

6.564E-01 

1 .3965 + 00 

2* 7686+0 1 

9.952E+03 

1.344E-01 

-1.180E-01 

— 6.395E-01 

1 .571E+00 

-2.8006*01 

9. 9395*03 

1v318E-01 

— 1 i433E-0 l 


1=9 

r- TH — P Q — PHE S t 

1 • I53E+00 0. 2.630E+01 1.003E+04 1.515E-01 1.117E-01 

- lVl 29E+00 ' 1.745F-01 2. 6675*01 — T.002E+04 "1.5055-0 1 8. 057E-02 

1.1045*00 3.491F-01 2.7045*01 1.000E+04 1.495E-01 6,OI5F : “02 

1 . 0555*00 6.981E-01 2.774E*0l 9.974E+03 1.475E-01 -8.690E-03 

r - 03 j pvOO 8.727E-01 • 2.808F* 01 9 .96 1 £ + 03 t .466E-01 - -3.7 1 1 E-02 

1.006F+00 1.0475*00 2.841F+91 9.948E*03 1.457E-01 -6.486E-02 

9”.' 8185- 0 1 17 222F + 00 2 78 7 35 *0 1 — 9.-9355*03 "1.4485-01 -9. 194E-02 

9.570E-01 1.396E*00 2.904E+01 9.923E+03 1.440E-01 -1.184E-01 

9. 3 25 E - 01 IV571E+00 2.935E + 01 9.91 IE *03 1.431E-01 -l;44l£“0* 


1 = 10 

— p r>r P 0 “PHE SI 

1.523E+00 0. 2.880E+01 9.985E*03 1.483E-01 9.8555-02 

r .4 9 1 E 4 ~00 1 ♦ 7 455- 0 1 27 ? O 3 5 * 0 1 : 9.9725*03 1 .4886-01 6r840t-'O2“ 

1 .4595*00 3.491 E. - 01 2.925F*01 9.9605*03 1.492E-01 3.969F-02 

I-4P6F * 0 0 5.2365-01 " 2." 9475*01 9. 9495+03 1 . 497E-0 1 1 . 1S1E-02 

1.3945*00 .6.9815-01 2.9695*01 9.937E+03 1.502E-OI -1.623E-02 

— 1.361 E ♦ 00 "8 .7 2 7 F. - 0 1 2.9905*01 ~ 9.926E+03 — 1 .507E-01 -4.333E-02 

1.3295*00 1.047E+00 3.010F+01 9.915E+03 1.511F-01 -6.979E-02 

1 .29 7 5*00 1 . 222ET 0 T) 3.030 6 *01 9.904E+03 1. S15F>01 -9.0G2E-02" 

1, 2645*00 1.3965*00 3.0505*01 9.894E+03 1.520E-01 -1.208E-01 

1.2326*00 1 .57 IE* 00 3.069F*01 9.883E + 0T — 1.5245-01 -1 .454E-01 
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H PHI PHO - GAM 

4.0486*07 8.0006-01 3.7716-06 1.1826*00 

4.05 7 E+ 0 7 8.000 r -0 1 3.8346-06 tTt826*i)fr 

4. 0676*07 8.000E-01 3.896E-06 1.1826*00 

4.0766*07 “ 8.0006-0r~379566-06 1.182E+00 

4.0856*07 8.000E-01 4.0156-06 1.181E*00 

4.0946*07 8.0006-01 4.0736-06 1.1816*00 

4.1026*07 8.0006-01 4.1296-06 1.181E+00 

4 ~ ;rr oE » (T7 — b twot^ p i — 4 n 84r-06 — rmri e ♦ oo 

4.1 19E*07 8.0006-01 *.238E-06 1.181E+00 

“4".126E*0 7 8.0006-01 — 4.290E-06 1.1616*00 


p PHI RHO GAM 

4.064E+07 8.0006-01 3.8466-06 1.182E*00 

“4.0756*07 8.0006-01 3.9096-06 1.1826*00 

4 • 086E*07 8.0006-01 3.971E-06 1.182E*00 

4. 0966*07“ 8.0006-01 4vG32E-06 1.18lE*00 

4 • 1066*07 8.0006-01 4.0926-06 1.1816*00 

4 • 1266*07 8.0006-01 4.2066-06 1.181E*00 

4.1366+07 8.000E-0 1 4.2616-06 1.181E + 00 

4.1456*07 8. OOOE-Ol 4.3156-06 1.180E+00 

4.1546*07 0.0006-01 4.368E-06 1.180E+00 


4 . 104E + 07 
4.1146*07 
4.1236*07 


PHI- 


8.000E-01 

8.0006-01 

8.000E-01 


-RHO 


4 • 046E-06 
4. 103E-06 
4. 1586-06 


- GAM 

1.181E+00 
1 • 18 1E + 00 
1 • 181E + 00 
1-.-1816+00- 


4.1416*07 8. OOOE-Ol 

~4~. 1496*0 7- 8.0006-01 

4 • 158E+07 8. 000E-01 

4.1666*07 8. 000E-01 

4.1746*07 8. OOOE-Ol 


4. 2656-06 1.1816*00 
4.3176-06 1 • 1 80E*00 
4.3676-06 1 • 1806*00 
4.4176-06 1 .1806*00 
4.465E-06 1.1806*00 


l-l-E-06— l-rie 0 E-+0 0 


H - PHI RHO GAM 

4.1616*07 8. OOOE-Ol 4.367E-06 1.180E*00 

^ . 1 -676*07 StOOGE-— 01 4-^4026—06 1-.-1 006*00 


4. 172E+07 
4.1786*07 
4.1836*07 
4.1896*07 
4.1946*07 


8.0006-01 
8.0006-01 
8.0006-01 
8.0006-01 
8. OOOE-Ol 


4.4366-06 

4.4706-06 

4.503E-06 

4.5356-06 

4.5666-06 


1.1806*00 
1 * 1 80E*00 
1.1806*00 
1.180E+00 
1 .1806*00 


4vlP 96-*fr? dvO-Q 0 E — 01 4-.-S 9 7 6-06 6.-180 6*- 0 0 


4.2046*07 8. OOOE-Ol 4.6266-06 1.1806*00 

4.2086*07 8. OOOE-Ol 4.6566-06 1.1796*00 
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R 

1 .922E+00 

— 1.08 IF* 0 0 

1 ,840E«-00 

1 .799F+00 

1 .759E*00 


11 

— -nr - 

o. 

1. 745E-0T' 
3.491E-01 
5. 2382-01 
6.981E-01 


3 , 204 F ♦Ol 
-3V204E+01 
3.205E+0 1 
“O'; 205F + 0 1~ 
3.206E+Q1 


0 

9.926E*03 
• 919E*0"3' 
9.910E+03 
,902E>03 
9.894E+03 
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PHf- 

1.553E-01 
1 .559E-01 
1 .565E-01 
i .570E-01 
1.576E-01 


-st- 


6.377E-02 
5.576F-02 
2.835F.-02 
1 .546E-03 
-2.4b6E-02 


H p+tf RH© GAH— 

4.236E+07 8.000E-01 4.782E-06 1.179E+00 
4.236E+07 8.000E-01 -4.78SE-06 1 . 179E+00 
4. 2362*07 8.000E-01 4.788E-06 1.179E*00 
4 ♦236E + 07 - 8.000E-01 - 4.790E-06 1.179E*00 
4.235E*07 8.000E-01 4.793E-06 1.179E*00 


1 .718E+00 
1 *677E*00 
1.636£*O0 
1 • 595E *00 
1 ,554E*00 


8.727E-01 
1.047E+00 
1T222E*00 
1 ,396E*00 
"T. 57 IE ♦0 0 


3.206F. + 01 
3.207E+01 
3.2Q7F+01 
3.207E*01 
3.208E*01 


9.8H6E+0 J 
9.878E+03 
9.870E+03 
9 . 663E *0 3 
9.856E+03 


l .bttit-Ul 
1 .586E-0 1 
1 .S91E-01 
1 .596E-01 
1 .601E-01 


-b.U^bh-02 
-7.S26E-02 
-9.966F-02 
-1.235E-01 
-1 .467E-01 


I = 


2.225E*00 
T. 1782*00 
2.130£*00 
2.083t*W 
2 . 036E *00 
i;989E*00 
1.94lE*00 
“ITS 9 AT*' O'O' 
1.847E*00 
TTT99EPnTTT 


12 
TH 

0. 

rv745r-ni 
3.491E-01 
5.236t:^rr 
6.981E-01 
8.727E-01' 
1.047E+.00 
-1.2222*0 0 
1 * 396E+00 
IT57T E>~mr 


3.462E*0l 9 

'3.443F + 01 9 

3.425F*0l 9 


'3.40 /F*01 
3.390E+01 
3T3T3F*trr 


3.357F+01 9 

“3.340F*0l 9 

3.325F+0 1 9 

-9 


‘3. JU9E + 01' 


,887E*03 

,88lE*03 

.875E+03 

,8b9t*03~ 

.863E+03 

7857 F *03 

,852E*03 

• 84 6 F. + 03 

.841E+03 

.836E*- 0 3 - 


— ph r 

1 .b56E-0 1 
1 .654F-01 
1.653E-01 
1 .852E‘ 0 r 
1 .650E-01 

l.b49E-0l 
1.648E-01 
1 .647E-01 
1 .646E-01 
1 . 645 F - 0~1 


— SI 
7.057E-02 
4.408E-02 
1.817E-02 
- 7 .17 f£- 03 


-3.195E-02 

--5.616E-02 

-7.980E-02 

-1.029E-01 

-1.254E-01 


4 . 2 Jbt*U ( 
4. 235£*07 
4 .235E+07 
4 .235£*07 
4.235E+07 


8.0UUt“Ul 

8.000E-01 

8.000E-01 

8.000E-01 

8.000E-01 


4.798E-06 
4.800E-06 
4 . 803E-Q6 
4.805E-06 


i • i fvt+uu 
1 • 1 79E + 00 
1 .179E+00 
1.179E+00 
1 • 1 79E + 00 


h — 

4.288E+07 
4.284E*07- 
4 • 280E+07 
-zrr2T6tr*07- 
4 • 272E + 0 7 
4.269E*07- 
4 • 265E + 0 7 
4.261E*07 
4.258E*07 


PHI 

8.000E-0I 

8.000E-01 

8.000E-01 

8.000E-01 

8.000E-01 

8.000E-01 

8.000E-01 

8.000E-01 


RHO GAM — 

5.106E-06 1.178E+00 

5.083E-06 1 • 176E*00 

5.060E-06 1 • 1 78E+00 

— 5vO37t:-06 lvl 78E ♦ 00 

5.015E-06 1.178E+00 

- 4.994E-06 1.179E+00 

4.973E-06 1 • 179E+00 

4.953E-06 1 . 179E+00 

4.933E-06 1 • 179E+00 


3E-06 — 1-.-179EKH)- 
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l 


I = 13 

— R TH :* > — "Q PHE S I 

2.225E+00 0. 3.4626+01 9.3886+03 1.712E-01 3.3366-02 

2 71 78E ^ ~0t ) 1 tT¥ 5E-01 — 3.4436+01 9 t P 4 ~ 3E *03 lv703E«01 1.1406-0 2 

2.1 30E + 00 3.491E-01 3.4256+01 9.496E+03 1.695E-01 -1.0096-02 

2 . 083E + 0 0 5.2366-01 3.4076+01 ” 9 . 54 66 *0 3 - "1 ; 6 8«F - 0 1 -3.1 1 IE-02 

2.036E+0O 6.981E-01 3.390F+01 9.599E+03 1.680E-01 -5.165E-02 

1 .989E *00 ' 8.727E-01"' 3.3736+01 9.649E+03 1 .673F-0r -7. 1 73E-02 


I = 13 


2.225E+00 
2vl 7 8 E + 0 0 
2 • 1 30E + 00 
2.0836+00 
2.036E+00 
n.9896+00 
1 .94 1 E * 0 0 
-178946 + 0 0 
1.847E+00 


— — TH 

0. 

lv745E 
3.491 E 


3.462E+01 
-01-3.44 3F*01 
-01 3 . 425E *01 

6.981E-0 1 •" 3.390E + 01 


Q- 

9.388E 

9.443E- 

9.496E 


+ 03 
♦ 03 — 
+ 03 
*03 


8.727E-01 
1 .047E + 00 


3 .3736+01 
3.357E+01 


1 .2226+00 — 3.-34 0E + O1- 
1 .3966 + 00 3 . 325E*0 1 


9.599E+03 
9.649E+03 
9.697E+03 
-9v745£ + 03 — 
9.791E+03 


— PHE 

1.712E-01 
1.-70 3E -01 
1 .695E-01 
tTt )88 E -0 - l 
1 .6806-01 
1 .6736-01 
1 .b65E-01 
lib58F-01 
1.6516-01 


— SI 
3.336E-02 
1.1406-02 
■1.009E-02 
>3.11 16-02“ 
■5. 1656-02 
■7. 1736-02 
-9.1 33E-02 
> 1 • 1 05E-0 1 
>1 .2916-01 


1.799E + 00 1.5716+00 

3.3096*01' 

~-9ir836t-*03 1;645E~ - Q 1" 

-1.4/JE-trt- 

I = 14 

P - TH — 

P 

Q PHE 

S 


2 .65 1 E* 00 
2.6o n^nrrr 
2.5506+00 
-2-.499E + 0O- 
2.4496+00 
-2v398E+D0 
2.3476+00 
-2T29-7E-rirO- 
2.2466+00 
2 rl 956 + 00- 


0. 

--T7T4sr-(r 

3.491E-01 

- 5.2366-01- 
6.9816-01 

8.7276-01 

1.0476+00 
— tri 

1.396E+00 

- 1.5716+00 


4.2446+01 

j | / *. j-r , A 1- 

^rrr99E *01 
4.1556+01 
4 .1 126 + 01“ 
4.0696*01 
4.0286*01 
3.9876+01 
-3r94-7E-MH“ 
3.9086+01 
-3 .870E-+01- 


9.370E+03 
-*^~3- 
9.4016+03 
-9.4166*03 
9.431E+03 
9.4456*03 
9.4596*03 
- 9 v 47 3£-+'03- 
9.4876*03 
-^isooe+'oa 


1.9506-01 
-hr" 
1.8856-01 
-1 .-8536-01 
1.8226-01 
1.791E-01 
1 .7626-01 
-M326-0+ 
1.7046-01 

-lv^76E-01- 


3 • 922E-02 
-3v(mE-fl2- 
2.2896-02 
I *4976-02 
7.2026-03 
-4.052E-04 
-7 .853E-03 
>+v5 146-02- 
•2.2276-02 
-2.9246-02 


I = 15 

— — R T H P 9 

3.4836*00 0. 6 . 0616+01 9.359E+03 

7- 3 .4266 + 00 1 .7456-01 5.9366*01 9*3 7 16^*03 

3.369E-00 3 . 4916-01 5.8126+01 9.384E+03 

-3.-311 E ♦ 0 0- 5.-2366-0 1 5.69 16 + 01 9.3976 + 03 

3.254E+00 6.9816-01 5.5726*01 9.4O9E+03 


r+4€ 54 

2.318E-01 4.217E-02 

-2.2846-01 3 .428E-02 

2.2496-01 2.6536-02 

2. 216E-01 - 1.8906-02 
2.1836-01 1.141E-02 
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3. 197E*00 

6.727E-01 

5.455E + 01 

9.421E+03 

2.150E-01 4.040E-0J 

3. 140E+00 

1.047E+00 

5.340F+01 

9.433F+03 

2.119F-01 -3.194E-03 

3.083E+00 

1 . 222 F +00 

5.227F+01 

— 9-.444E + 03 ‘ 

2.0A7E-01 -1.030E-02 

3.025E+00 

1 .396E+00 

5.117F+01 

9.455E+03 

2.057E-01 - 1 . 727F.-02 

2.968E+00 

1 .5 7 1 E ♦ 0 0 5 . 0 0 8 F ♦ 0 1 

9.467E+0 J 

2 . 026E-01 -2.411F-02 

I 

= 16 




R 

TH 

P 

Q 

PHE si 

3.483E+00 

0 . 

1 .000F + 01 

9 . 900E ♦ 03 

4.650E-14 -5.06SE-15 

3 .426E+00 

"i; 745E-0 1 

’ 1.477F+01 

9.843E+03 

2.41 IE-02 -2.869E-03 

3.369E+00 

3.491E-01 

1 .946F+01 

9.798E+03 

4.783E-02 -5.691F-03 

3.“3liE>00 

^5T23'6e-0T“ 

”2. '40 7F' + 0l'' 

9. 748E + 0 3 

f»ll4t-02 -8.4&feE-03 

3.254E+00 

6.981E-01 

2.860F+01 

9.699E+03 

9.406E-02 -1.1 19E-02 

3.197EV00 

" 8V727F-01 

3.306F+01 

9 ,651 E+ 03 

1 . 166E-01 -T.387E-02 

3.140E+00 

1 .047E+00 

3.743F+01 

9.603E+03 

1.387E-01 -1.650E-02 

3.083E + 00’ 

" '"IV222E + 00 

"4.173F>or 

9.557F + 03 

1 .604E-01 -1.909E-02 

3 • 025E ♦ 00 

1 • 396E+00 

4.594F+01 

9.51 1E + 03 

1.817E-01 -2.163E-02 

2T968E + 0 D” 

1 .57 IE + 00 

» UOSE + Ol 

9 . 467 E ♦ 0 3 

2 •026E-0 1 -2*4 r 1 E-02 

I 

= 17 




R 

““ m 

p 

0 

_ F>H£ 5 j - 

2.352E+01 

0. 

1 .000F + 01 

9 • 90 0E + 03 

0. 0. 


77T6 T ?ETTn r7755T=in 1 .000F+01 FTWOF+TO a . 7 3 2F - 2 9~- -2Tl 6 0 F-ltr 

2.385E+01 3.491E-01 1.000F+01 9.900E+03 1.753E-28 -4.335E-16 


“ 2.401E + 0I “ 5.236F-01 ‘ T.O00F*01 * 9V90 0E+03 “‘2S638E-28- -6.525E-16 
2.41PE+01 6.981E-01 l.OO'OE + Ol 9.900E»03 3.529E-28 -8.730E-16 

— 2 .434E + 0 1 " 8 .727E-01 1.000F+01 9.900E+03 4.427F-28 -1.095F-15 

2.451E+01 1.047E*00 1.000F + 01 9.9Q0E+03 5.3 31F-2 8 -1.319F-15 

2 ^484E + 0 1 1 . 396E+00 1.000E+01 9.900E + 03 7.1S6E-28 -1.770E-15 

2 *'50 0 E + 0 V 1 .57 1 E ♦ 0 0 — 1 0 0 0F+ 0 1 9.9 0 0 E+3)3 8 ; 07 8 E - 28 - 1.99 8 E - 1 5 
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1.4536+Uf 0. 1.9bbt-Ub 

1. 4436*07 0. 1 . 943E-05 

1.432E+07 0.- 1.921E-05 

1.422E+07 0. 1.900E-05 

I . 4 1 2 E ♦ 0 T 0 • 1. 879E-05 


1 . J4Vt>UU 

1.350E+00 

1.351E+00 

1.351F+00 

1.352E+00 


PHI — RHO - GAM 


9.950E+06 

1.045E+07 

1.0936+07 

0. 

~0. 

0. 

6 *456E-06 

7.924E-06 

9.368E-06 

1.387E+00 

1.383E+00 

1.379E+00 

1.141 E+OT 

0 * 

' rrOT9E=95— 

~t » 3756^+00 

1 .1896 + 07 

0. 

1.218E-05 

1.371E+00 

- 1-.235E + 07 

0. 

1; 3556-0 5 

1.367E+00 

1.2806+07 

0. 

1.490E-05 

1.363E+00 

1.3256+07 

— 0. 

1-V622E-05 

1.3596+00 

1.3696+07 

0. 

1.752E-05 

1 . 356E+00 

1 •412E+07 

0 • 

fv8~79E-05 

~tr35-2E+0fr 


H PH I RHO GAM 

9.950E+06 0. 6.4566-06 1.387E+00 

9*9506 + 06 Ov 6v45y6tr^6 lrv38 7E*GO 

9*9506+06 0. 6.456E-06 1.387E+00 

• 9 .950E+G6- 0 * ~6 .456E-06 1 . 387E + 00 


9.950E+06 0. 6.456E-06 1.387E+00 

- 9.9506+06 0. 6*4S6E-06 1.387E + 00 

9.950E+06 0. 6.456E-06 1.387E+00 

~~~9 . 95 ~9E~^96 0 - ; 6t 456E-06 — t~.-3 87E+00 

9«950E+06 0. 6.456E-06 1.387E+00 

-9T950E + 06 - 0^ 6*4566-06 1.387E+00 
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KOUM a 0 


X 9 

2.749Q0E+01 

- . . .. . — — 

TABLE Ilia 





J * 1 



THETA B 

0* 

I 

Q 

M 

Q 

PHE 

si 

1 

1 •5C0t+vl 

3, (Hi OF *02 

7,*?OE+03 

1.745E-01 

0. 

a 

5 .510E+01 

r. f-UOF + 02 

7.420E+03 

1.745E-01 

0. 

3 

1.520E+Y1 

3 . . .. 0 o P + 0 2 

7.420E+03 

1.745E-01 

0, 

4 

1 -530t+yl 

3 . 0 0 0 F + 0 1 

7.420E+03 

1 .745E-01 

0. 

3 

1 •5&Ok+0 1 

3 # ()o!>r + 02 

7.420E+03 

1.745E-01 

0. 

b 

1 •550fc.+ul 

3.00 Of +02 

7 *420E+03 

1.745E-01 

0. 

7 

1.56nb*0l 

n . 0 (> U 3 + U 2 

7.420E+03 

1.745E-01 

0. 

JM 

1 • 57 Ofc. + 01 

F.f.uor: +02 

7.420E+03 

1.745E-01 

0. 

9 

l + o 1 

S.OOOh *02 

7.420E+03 

1.745E-01 

0. 

10 

1 »5«>0t + wl 

5*(iUOt + 02 

7, 4?0E+03 

1.745E-01 

0, 

11 

1 *600fc+01 

3. OOOr. +02 

7.420E+03 

1.745E-01 

0. 
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M 

H 

PHI 

2.211*>00 

6.450E*07 

e.ooop-oi 

2.211F+00 

6,A5UE->07 

6.000F-01 

2.2UE + 00 

6,450E*07 

8,000F-01 

2.2llP*00 

6«45Q£*07 

8.000P-01 

2.2UP + 00 

fe.45U£*07 

8.000F-01 

2 • 21 1 F + 00 

6,45OE*07 

0.OOOP-O1 

2.2MP + 00 

6.450E*07 

8,oooe-oi 

2.211P+00 

6.450E+07 

8.000P-01 

2#21i>+00 

6 • 45 ♦ 07 

e.ooot-oi 

2.211P+00 

6.450fc*07 

e.oooF-oi 

2.211P*00 

6,450t*07 

e.OOOE-Ol 
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HHO 

OAM 

T 

5.089E-05 

1 • 147E*00 

4.728E+03 

5.089E-05 

1 . 147E*00 

4»728£*03 

5.089E-05 

1 * 1 47E + 00 

4.728E+03 

5.089E-05 

1 . 147E+00 

4*728E+03 

5.089E-05 

1 ,147E*00 

4.728E+03 

tj,089E-05 

1 « 147E+00 

4.720E+O3 

5.089E-05 

1.147E+00 

4.728E+03 

5.089E-05 

1 .147E+00 

4.720E*O3 

5.089E-05 

1.147E+00 

4.728E+03 

5.089E-05 

1.147E+00 

4.728E+03 

5.089E-05 

1 ,147E*00 

4,728E*03 



166. r-Yol . P*ge 164 

KOUNT = 110 X ■ 3.07600E+01 

. TABLE _LLLb 


j.« i 

I R P 

1 1*454E*«1 1 • 77 OC*02 

-? I.482E+01 1. 700F*02 

3 1.507E*wl 1.640F+02 

4 1 *532t*0l l.*0*F*02 

5 -1.555fc*wl 1 * 5 s>7f; + 02 

6 l*b78t*ol l.bl3E*02 

7 l •600t*wl 1.63«F*02 

■8 1.6?lE*wl l.f>b«E*Q2 

9 l.b4?t*vl 1.66bF*02 

10 l,663t+wl 1.8S7 e*02 

1-1 l*685fc+Ul 1.6318*02 - 


THETA ■ 0« i — 

Q PHE SI 

3.611E+03 -1.551E-01 0. 

8,650E*03 -1.149E-01 0. 

b , 686E+03 -b.491E-02 0. 

b.707E*03 -b,567E-02 0, 

6*71 1E+03 -2.279E-02 0, 

b« 702E + 03 1.493E-02 0, 

b # 687E + 03 5.729E-02 0. 

3.675E+03 1.037E-01 0. 

b # 67lE*03 1.&36E-.01 0. 

8.676E+03 2,064£-0i 0. 

3.691E+03 2.618E-01 0, 


. 

J 5 2 

- 

‘ 

-THETA a 

8.7P600E-02 

I 

R 

P 

Q 

PHE 

si 

1 

1.46fet*01 

1.74&E+02 

d,586E*03 

-1.604E-01 

-1.177E-01 

-2 - 

1 *49 It + 0 1 

1 *705F*02 

8.627E+03 

-1.330E-01 

-7.593E-02 

3 

l.5l4b*^l 

1 ,68br + 02 

b.S40E+03 

-1.020E-01 

-7.00BE-0? 

4 

l,53*>t + Ul 

1.709F+02 

b.62SE+03 

-b .293E-02 

-7, 150E-02 

% - 

1 *5$7t*l* 1 

1.777F+02 

tf *58fcE+03 

-1.800E-02 

-7,pS2E-02 

6 

l*577t+ol 

1 ,67bE>02 

B.532E+03 

2.862E-02 

“7.073E-02 

7 

1 «59t)t + w 1 

1.9898*02 

8.474E+03 

7.452E-02 

-6.6J7E-02 

6 

l*4J*t*ul 

2.0948*02 

6.424E+03 

1.197E-01 

-6.096E-02 

9 

l*h3lt+ol 

?. 170F + 02 

8.389E+03 

1.653E-01 

“5* B2OE»0? 

10 

l,64*E*0l 

2.2038*02 

8.373E+03 

2.123E-01 

“5 , 747E-02 

n 

1.66St*ol 

2.1898+02 

d,37BE*03 

2.61BE-01 

-5.932E-02 


t . “J 



THETA a 

1.74b20E-01 


sj ^ «? 



i 

R 

P 

Q 

PHE 

SI 

l 

1.479t*01 

1 . 7obE+02 

8.609E+03 

-1.656E-01 

-1.113E-01 

2 

1.502t*ol 

1 .744F + 02 

8.607E+03 

-1.416E-01 

-7.177E-02 

3 


1 . 8558*02 

d,548E*03 

-1.012E-01 

-6.276E-02 

4 

l*54lt*«Jl 

2-0148+02 

d • 466E ♦ 03 

-4.886E-02 

•b«90bt ,, 0? 

5 

1 *S 5 *t +01 

2 . lobF*U 2 

8.383E+03 

7.782E-03 

-5.563E-02 

6 

1 • 573 t + v ' 1 

? » SS't’F * 02 

8.306E+03 

6.295E-02 

-5.146E-02 

7 

1 « 5P9t * o 1 

2 , 52 ir *02 

b« 234£+ 03 

1.129E-01 

-4.632E-02 


1 ,ftO?t+vl 

?.*'bb c *02 

8, 166E+03 

1.559E-01 

-4 # o«3E-02 

9 

l*6l?t*vl 

? . b 33 K + 02 

d, 107E+03 

1.930E-01 

-3.420E-02 

10 

1 . 63 1 1 * ^ 1 

?. 9 iF 8+02 

o, 068E+03 

2.269E-01 

-3.B41E-02 

il 

1 ,645t»0 1 

2 . 9501 : *02 

8.061E+03 

2.&18E-01 

-3.973E-02 


il 

j 


„ 

THETA * 

2.61780E-01 

I 

P 

M 

Q 

PHE 

si 

1 

1 * 49 1 1 ♦ 0 1 

1 . 7 U 3 F +02 

b.bOlE+03 

-1.709E-01 

-1.217E-01 

2 

1 , 5 l 2 t + 0l 

2.0318*02 

d,458E*03 

-1.117E-01 

-6.225E-02 

3 

1 * 52 Ht*u 1 

?.* 9 Br*u 2 

b.246E*03 

-4.384E-02 

“4. J 09E-02 

4 

1 * 542 t + 0l 

P.^bF + O* 

d, 049E + 03 

1.967E-02 

“3.086E-0? 

S 

l,555E*ol 

3.4 10C- + 02 

7.B96E+03 

7.456E-02 

-2.5J1E-02 

6 

1 *567t * 0 1 

3. 700**02 

7.797E+03 

1.227E-01 

-2.197E-02 

7 

1 *57Bt>0l 

3 . hbJE + 02 

7.750E*03 

1.6S6E-01 

-1 .997E-0? 

% 

1.5fi9t+ul 

3, ht>oF + U2 

7.746E*03 

2.031E-01 

-1 *90«E«0? 

9 

1 . 60 lE*wl 

3 ,mu 5 F *02 

7.746E+03 

2.332E-01 

-1 .Q2S.E-0? 

10 

l.<S13 t *ol 

3.74Qr+C2 

7 , 786E + 03 

2.528E-01 

-2.2S3E-02 

l ; 


~ 7 ,, f 4 0 - 

7 VOuf 4 r> "i 

. t. 1 m ZT — * 1 

- - *’ H oS - ^ -<• 
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M 

H 

2,734F*00 

5.496E*o7 

2,753f*00 

5.46lE*07 

2.771F*00 

5«431t.07 

2, 76?f*00 

5.41<;e*07 

2.784F+00 

5.40oE*07 

2,779F*00 

5.41 /E+0? 

2.772F+00 

5.430E*07 

2.766F*00 

5.4*0£*07 

2.763F *00 

5,444£ + 07 

2,766r*00 

5,43y£+o7 

2.774F*00 

5.42bE+07 


M 

H 

2»743f+00 

5.467£ + o7 

2.752F+00 

5.46uE*07 

2.75Bf*P0 

5.452E+07 

2.75lF*00 

5.46bE*07 

2.73?F*00 

5»49«t*o7 

2.706F*00 

5.54SE+07 

2 »677f*00 

5.596E*07 

2.652f+00 

5.642E*o7 

2 » 634F *00 

5.67JE+07 

2*627F *00 

5.666E+07 

2.63OF.00 

5.66 1E+07 


M 

H 

2 • 75?F+00 

5.449E+07 

2.741f+00 

5.479E+07 

2.711P+00 

5.534E+.07 

2»67 jF* 0 0 

5.fcoe>t*o7 

2.63 IF* 00 

5 . 676E+07 

2.594F+00 

5.745E+07 

2.560 r *00 

5.806£+07 

2.529c + 00 

5.863E+07 

2.502F+00 

5.9 12£*o7 

2.4B5F+00 

5.943E«-07 

2,482r 4 00 

5.949E*n7 


M 

H 

2.75^f*00 

5,44bE*07 

2,667f + o0 

5.612E+07 

2.565F+00 

5.79 /E + f)7 

2.475F+00 

5.961E+07 

2.409F400 

6 , 06 JE*07 

2. 366f*00 

6. 162E+07 

2. 346F4QQ 

6. 190t*(*7 

2. 344C.00 

6,c02fc.07 

2. 353r *00 

6. lfce£*07 

2. 36le-*00 

b. 1 7 0 £*o7 


a n a «_ i u 


PH! 

RHO 

8.000E-01 

2.085E-05 

8.000F-01 

1.995E-05 

8.000E-01 

1.934E-05 

8.000E-Q1 

1.090E-O5 

6.000F-01 

1.890E-05 

8 . OOOE-0 1 

1 , 907E-05 

8.000F-01 

1.932E-05 

b.oooe-oi 

1.952E-05 

8.000F-01 

1 .960E-05 

8.000E-01 

1.951E-05 

8.000E-01 

1.925E-05 

PHI 

RHO 

S.OOOe-OI 

2.03SE-05 

8.000E-01 

1.999E-05 

S.OOOfr-Ol 

1.980E-05 

8.000F-01 

2.004E-05 

6.QQ0E-01 

2.072E-05 

8.000E-01 

2.172E-05 

8.000E-01 

2.285E-05 

8.000F-01 

2.389E-0S 

8.000E-01 

2.464E-05 

e.oooF-oi 

2.496E-05 

8.000E-01 

2.483E-05 


PHI 

Rho 

b.oooe-oi 

1.995E-05 

8.000E-01 

2.O30E-O5 

8.000E-01 

2.1&1E-05 

8.000E-01 

2.309E-05 

8.00QF-01 

2.479E-05 

e.oooF-oi 

2.644E-05 

8, OOOf-01 

2.0O5E-O5 

8.0GOF-01 

2.963E-05 

8.000E-01 

3. 105E-05 

B.OOOc-01 

3.199E-05 

8,oooe-oi 

3.215E-05 


PHI 

RHO 

8.0QOF-01 

1.992E-05 

8.000F-01 

2.325E-05 

8.000F-01 

2.783E-05 

8. 000E-01 

3. 251 E- 05 

8.000F-01 

3.647E-05 

8.000F-01 

3.921E-05 

8.000E-01 

4.056E-05 

8.000P-01 

4.069E-05 

8 . OOOF-0 1 

4.012E-05 

3 . 000p“0 1 

3.952E-05 

U r> m or « m t 

O *4 ? ^ ^ ^ 


6AM 

T 

1.5158E + 00 

4.248E*o3 

1.158E400 

4.?31E*03 

1,1S9F*00 

4.216E+03 

1.159E+00 

4.207E+03 

1.1696*00 

4.205E+03 

1. 159E+00 

4.209E.03 

1 . 159F+00 

4.216E+03 

1.159E+00 

4.221E*63 

1.159E*00 

4.223E+03 

1.159E*00 

4.220E+03 

1.159E+00 

4.214E403 


GAM 

T 

1.158E+00 

4,?36E*03 

1. 158E+00 

4,231E*03 

1 ♦ 158E+00 

4.227E*03 

1 , 1 58E+00 

4.P33E.03 

1 . 158E+00 

4.?49E*03 

1,1 57E+00 

4.272E+03 

1 , 156E + 00 

4.297E+03 

1.156E400 

4.319E*03 

1 • 155E*00 

4.335E*03 

1 , 155E+00 

4 # 34 1 E*03 

1.155E+00 

4 . 338E+Q3 


GAM 

T 

1.159E+00 

4.227E*o3 

1.158E+00 

4.240E*03 

1 . 157E+00 

4.267E+03 

1. 156E*00 

4.302E*03 

1 , 155E+00 

4 * 337E*03 

1.154E*00 

4.370E+03 

1.154E+00 

4.401E«03 

1 . 153E *00 

4.429E*03 

1, 152E+00 

4.454E+Q3 

1.152000 

4.469E*03 

1 , 152E+00 

4.472E+03 


GAM 

T 

1 , 159E*00 

4.226E.03 

1 ( 1S6E«00 

4.305E.03 

1, 1 54E+00 

4.396E+03 

1,152E*00 

4.478E*03 

1 . 150E+00 

4.540E403 

1.149E400 

4.580E+03 

1, 149E*00 

4.599E403 

1.149E*00 

4.600E.O3 

1.1495*00 

4.592E+03 

1,).49e*00 

A,584e*03 

1 4 /, 

4 . CJ 3 4 ^ < 4 ^ 
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I 

J = f> 
P 

p 

Q 

THETA st 
PHE 

3,490406 

SI 

1 

l.?02t*^l 

2 . bl OF + 02 

8.179E+03 

-1 .9716-02 

0. 

? 

l«Sl4£«ul 

■*,41^f + 06 

7.733E+03 

9,7716-02 

0. 

3 

1 »5? ( ?t + vl 

4 • 6736 + Qd 

7,s37E*03 

1 • 484E-0 1 

0. 

4 

1*5356*01 

A.S< ,96*02 

7.460E+03 

1.677E-01 

0. 

5 

1 *5456+^1 

4,92fcr+02 

7.4396*03 

1,7286-01 

0. 

6 

1 ,5556+01 

4,44bF+06 

7.434E+03 

1.738E-01 

0. 

7 

1 *56' Pt + O 1 

4,9506 + 1)2 

7 .4336*03 

1.741E-01 

0. 

ft 

l»575t*ol 

4.9*-<?F. + 0£ 

7.4366*03 

1.7506-01 

0 , 

9 

1.5P5K+01 

4.6 e 2r+02 

7.452E*03 

1.795E-01 

0 . 

10 

1 *595t *0 1 

4,69(.'F + 02 

7.504E+03 

1.935E-01 

0. 

11 

1 »6o5t*w 1 

4, ji*9f+02 

7.6136+03 

2.225E-01 

0. 
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<7U' • i j 


M H 

2.53lr*00 5.8By£+07 

2.338r*00 6.21JE+07 

2.257E+00 6.36oE*o7 

2.227F+00 6.420E*07 

2,2l8r*00 6.43b£+07 

2.216r*00 6 . 43*rf£ + o7 

2.216F + 00 6 » 44u t ♦ 07 

2.217F+00 6.43«E*o7 

2.223F+00 6.a?7E*07 

2.244F+00 6.387£ + o7 

2.289F+00 6.30SE*07 


PHI RHO 

S.OOOfr-oi 2,9^«F-05 
e.OOOF-Ol 4, H0E-05 
8.000E-fil 4.708E-05 

8,000f-01 4.955E-05 

8.000F-01 5.025E-05 

6.000F-01 5.041E-05 

8.000F-01 5.045E-05 

8, 00 0E*0 1 5.037E-05 

8.000F-01 4.984E-05 

8 1 OOOE“0 1 4.813E-05 

8,000F-01 4.470E-05 


GAM T 

1.153f*00 4,426E*03 

1.149E*00 4 ,606E*q3 

1.147E+00 4.683E*o3 

1.1A7E*00 4.712E+03 

1.147E400 4.721E*03 

1 « 1 47E + 0 0 4,722E*03 

1 • 147E *0 0 4.723E + 03 

1.147E*00 4.722E*03 

1 • 147E+00 4.716E*03 

1.147E+00 4,696E*03 

1 • 148E+00 4 • 653E*03 
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V. CONCLUSIONS 

A numerical method based on reference plane characteristics 
has been developed for the calculation of steady three dimen- 
sional supersonic nozzle-exhaust flow fields. This reference 
plane technique by incorporating three separate coordinate 
systems, permits the calculation of geometrically complex flow 
fields using a relatively small amount of computer time. 

The program developed has second order accuracy. The scheme 
used provides for convergence of reference plane properties prior 
to updating the cross derivatives. This is consistent with stan- 
dard practice in computational schemes employing characteristic 
grids. Further improvements in computational accuracy are 
achieved through the use of integral correction factors for mass 
flow and total energy. 

For end modules the external corner expansion is treated by a 
power series technique as described in this report. A further 
improvement in the program would be accomplished by replacing 
this expansion with a locally cylindrical reference plane sys- 
tem centered about the corner point. 

The program developed computes the entire flow field given ap- 
propriate initial data, external flow data and the vehicle 
geometry. All required coordinate transformati ons are perform- 
ed automatically. Shocks, contact surfaces and streamline pro- 
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jections are traced in individual reference planes to facilitate 
machine plotting. Additionally, initial data may be prescribed 
along a sweptback surface for a cylindrical or Cartesian co- 
ordinate system. Shock wave formation due to coalescence is 
detected, although these additional shock waves are not carried 
along as discontinuity surfaces. 

The program developed can analyze a wider variety of problems 
than described in this report. This geometric flexibility is 
provided by incorporating three (3) separate coordinate reference 
systems and tha ability to compute discrete surfaces of disconti- 
nuity in these systems. A program of this size and complexity 
requires a thoughtful selection of the appropriate reference 
systems, initial and external flow data, and body geometry poly- 
nomials. However, it Is felt that once the present program has 
been run and understood, modifications can be performed to yield 
program options encompassing a wide array of three dimensional 
supersonic flow problems. 
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APPENDIX 

CURVE FITS FOR r, h and p 

The variation of r (the equilibrium value of y) as a func- 
tion of temperature (T)» pressure (P) and equivalence ratio 
($>) is presented graphically in Figures (A1 ) , (A2 ) and (A3) 
from values tabulated in Reference (4). In Figure (A1 ) it 
can be seen that r is a strong function of T over the tem- 
perature range of interest, while the effect of varying com- 
position is small by comparison. Moreover, Figure (A2 ) in- 
dicates that r is moderately sensitive to pressure and the 
degree of sensitivity increases substantially as the tempera- 
ture level increases and dissociation effects become impor- 
tant . 

As a result of these observations, temperature is the pri- 
mary independent variable, while pressure is the secondary 
independent variable and composition acts as a perturbation 
variable. Thus, we can fit the function r ( T , P ,$ ) with a 
polynomial in T and add on a temperature dependent correc- 
tion term for the effect of pressure and a temperature in- 
dependent correction term for the effect of 

* 

An examination of Figure (Al) suggest that the function r (T) 
can best be curve fit by breaking up the temperature range 
into three intervals such that the function can be represent- 







PRESSURE (PASCAL) 


FIGURE A2 . r VARIATION WITH PRESSURE 
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ed by a parabola in each range. Choosing p = 10 pascal and 
<j>=l as our base, we therefore find three functions 


r a (T , 10 5 , 1) = - 1.833xl0“7 t 2 + 7.5X10’ 5 T + 1.367 
r 2 (T,10 5 ,l) = 2.0X10- 8 T 2 - 1.38xl0' 4 T + 1.423 
r 3 (T,10 5 ,l) = 7 . 27xl0" 8 T 2 - 4.57X10- 4 T + 1.85 


and define the basic temperature function as 


r(T,10 5 ,1) 


rid.io 5 ,!) 

- / r 2 (T,10 5 > l) y for 

r 3 (T ,io 5 ,1) 

\ J / 


T < 5 0 0 0 K ) 
< 500<T£2000°K 
T>2000°K 



Figure (A3) indicates that 3T. is constant in the two ranges 

a$ 

$<1 and $>1, but is a function of T. Fitting the function 
ar in each of the ranges of <j> we obtain 


ar 

a$ 


"nj (T)" 


r 

$<i 

for 


. n 2< T L 




where 

n : (T ) = 4xl0“ 9 T 2 - 2xl0‘ 5 T - 0.019 

n 2 (T) = 3.39xl0' 2 T°‘ 5 - 3.9lxl0 -4 T - 0.681 
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This now defines r as a function of both temperature and $ 
by means of the equation 


r(T,io 5 ,$) = r(T,io 5 ,i) + ($-1) AL ( 8 ) 

3 $ 

Finally, the effect of pressure must be included. From 
Figure {18} we observe that r may be approximated as 

r(T,P,$) = r(T,10 5 ,*) + m [1 og 10 (pxl0 5 ) -5] (9) 

where m is a function of T. Deriving m, we find 


m = 


■2. 15x10“ 8 T^+0.91xl0 _4 T- 0.0695 



1 T<1000°K 

> for < 

T>1000°K 


1 


(10) 


Summarizing, the final function obtained is 


r(T,P,$) = r(T,10 5 ,1) + m(l£-| -5) + !! ($-1) 


( 11 ) 


where the functions r(T, 10^,1), 9T_ and m are given by Equa- 
tions (4), (5) and (10) respectively. 

The curve fit for enthalpy is derived in a similar way. Figures 

% 

(A-4) and ( A - 5 ) present the variation of h with temperature, 
pressure and equivalence ratio. As was the case for r, the 




T(°K) 


FIGURE A4 . ENTHALPY AS A FUNCTION OF TEMPERATURE ( p=10 5 pa . ) 
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PRESSURE (pa) 


FIGURE A5. ENTHALPY AS A FUNCTION OF PRESSURE. 
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function h(T, <j>, p) is fit by a quadratic function of T, the 
coefficients of which are functions of <j> and an additive term 
for the effects of pressure. The resulting curve fit is sum- 
marized below. 

I h ( T > » 1 0 5 ) 

h(T,4>,p) = S 

h(T,4> ,p) 

where h(T,4>,p) - h(T,(p,10 5 ) 1 + 


(13) 

The basic function h(T, 4 >, 10 ^} is defined as 


5 > 2 - 275 ( 


ln p 

2.3 


■51 


for 


T^2000°K 


T>2000 K 


( 12 ) 


h ( T ,d> , 10 5 ) = 10 6 (a,T 2 +b 1 T + c 1 ) 


(14) 


with the coefficients aj, bj and c ^ defined below: 

For T < 2000°K and <J> < 1 

aj = 10" 7 ( - . 1 042cj> 2 + .82424) +.987) 

b x = 10“ 3 (.01 167cp 2 + . 15034) + .938) (15) 

Cj = - . 0284c> 2 + > .6731 4> + .4293 
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For 

T < 

2000°K and M 1 



a l 

= 10” 7 ( 1.787(f) 2 - 5.48(f) + 5.4) 



b l 

= lO -3 {-.1867(f) 2 + 1. 11$ + .176) 

(16) 


C 1 

= - .0933(f) 2 + 3.975(f) - 2.808 


For 

T > 

2000°K and <j> < 1 



a l 

= 10“ 6 ( 1 . 792(f) 2 + .3983(f) + .310) 



b l 

= 10" 3 ( -9 . 05(f) 2 - .07917(f) + .245) 

(17) 


C 1 

= 10 . 86(f) 2 - . 1 183(f) + .970 


For 

T > 

2000°K and <J> > 1 



a l 

= 10" 6 (4.81(f) 2 - 13.9(f) + 11.59) 



b l 

= 10‘ 3 (-23.08(f) 2 + 66.82(f) - 52.61) 

(18) 


Cl 

= 27.05(f) 2 - 73.73(f) + 58.39 



When the inverse function T(h,<J>,p) is required, it is obtained 
by an iterative solution of Equations (12) through (18) 

¥ 

The density is found by obtaining a curve fit for the mixture 
molecular weight and using the equation of state 
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P = M (19) 

RT 

where R is the universal gas constant and m is the molecular 
weight. 


The behavior of m with T,p and <£ is illustrated in Figures (A6) 
and (A7). We see that for temperatures less than 2000°K» m is 
essentially independent of temperature. The discontinuity- in 
slope of m( ) shown in Figure (A 7) requires that the equivalence 
ratio range be split in two. Thus, 

for T s 2000°K 




1 . 53<p 2 -5 ,895(J) + 28 .965 


i — i 

V 1 
-0- 

< 

* f or < 

1 .6O<j)Mo,60 + 33,6 


$ > 1 


( 20 ) 


For the higher temperature range, it is convenient to employ 
the form 


m = mU) - <s(P ,<t> ,T ) 


( 21 ) 


where 


<5 = d 2 ( P » $ ) 


( T-200o/ 2 ^ 
V ' 1000 ' 


( 22 ) 




70 


o> 

cr> 



FIGURE A6 , MOLECULAR WEIGHT AS A FUNCTION OF TEMPERATURE AND PRESSURE . 


Page A12 








FIGURE A 7 . MOLECULAR WEIGHT AS A FUNCTION OF EQUIVALENCE RATIO FOR T<2000°K, 
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and 

d 2 = a 2 (^ f ) 1 * 5 + b 2 (|^) + c 2 (23) 

= -2 . 3(£> 2 + 4.014? + 1.736 

b,. = 8 . 6 1 <}> 2 - 15 . 4 2 4? - 6.66 

For 0s<Ml * (24) 

c ^ = - 16 . 884? ^ + 33.21<j> + 14.53 

n 2 = . 4 3 7 5 cj> ^ + . 06254? + 2.08 

a 2 = -.322<f> 2 + 2.3634? + 1.905 
b 2 = 2 . 7 64> 2 - 7.56(f) - 3.68 

and for 1 sd? < 2 ? (25) 

= 3 . 6<J) + 7.364) + 27.15 

n 2 = . 47(f) 2 + 1 .825(f) + .350 



